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Abstract. We describe new computer-based search strategies for ex¬ 
treme functions for the Gomory-Johnson infinite group problem. They 
lead to the discovery of new extreme functions, whose existence settles 
several open questions. 
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1. Introduction and statement of results 

1.1. Group relaxations and extreme functions. Cutting planes are 
widely used in the state-of-the-art integer programming solvers. Impor¬ 
tant sources of general-purpose cutting planes are the master finite group 
relaxation of an integer program, which was introduced by Gomory in 1969 
[18] . and the infinite group relaxation by Gomory and Johnson [191120] , Due 
to the pressing need for effective cutting planes, the group problem has re¬ 
ceived renewed attention in the recent years, since it may be the key to new 
multi-row cutting plane approaches that have better performance than the 
ones in use today. 

Computer-based search has been used for the investigations of Gomory’s 
group problem and Gomory-Johnson’s infinite group problem since the very 
beginning, leading to the discovery of many cutting planes. In this paper, we 
develop new computer-based search strategies to carry forward the discovery. 

We restrict ourselves to the single-row (or, one-dimensional) problem. 
That is, we focus on only one row of a simplex tableau of an integer program. 
Suppose the row corresponding to some basic variable x is of the form 


x = ~ f + 

j=1 ( 1 ) 

x E Z_|_, 

Vj e ^+,Vj E {1,2,... ,m}, 


where denote the nonbasic variables. We assume / £l\Z, that is, 

the basic variable x is currently fractional. 

When all data is rational, there exists some integer q > 0 such that rj E 
X 7L for any j E {1, 2, ..., m) and / E ^Z. Gomory’s master finite (cyclic) 
group problem of order q is obtained by relaxing the basic variable x E Z + 
to x E Z and by introducing variables y(r) E Z + for every r E ^Z. Using 
the quotient group G/ Z, i.e., reducing modulo 1, and standard aggregation 
of variables whose coefficients are the same modulo 1 (see [9j Remark 2.1]), 
the relaxation of 0 takes the form 

Y r y( r ) = / + z > 

reG/Z (2) 

y(r) E Z + ,Vr E G/Z, 


where G = and / is a given element of G \ Z. The master finite group 
problem only depends on the parameters / and q, but not on any other 
problem data. 

Gomory-Johnson’s infinite group problem is obtained by further intro¬ 
ducing infinitely many new variables y{r) E Z + for every r E M. Formally, 
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again by aggregation of variables, it can be written as 

Y r y( r ) = /+ z > 

reG/z (3) 

y: G/Z —> Z + is a function of finite support, 

where G = M and / is a given element of G\ 7L. The infinite group problem 
only depends on the parameter /. 

We study the convex hull Rf(G / Z) of the set of functions y: G/T, —> Z + 
satisfying the constraints in ([2]) and in ^ for the finite and infinite group 
problems respectively^ The elements of the convex hull are understood as 
functions y: G/ Z —>■ M_|_. 

A function n: G/Z —)• R is called a valid function for Rf(G/ Z) if 

Y 7 T ( r )y( r ) ^ 1 ( 4 ) 

reG/Z 

holds for any y E Rf(G/ Z). Minimal (valid) functions are those valid func¬ 
tions that are pointwise minimal. Let IIf(G/Z) denote the set of minimal 
functions for Rf(G/ Z). Extreme functions are those valid functions that are 
not a proper convex combination of other valid functions. We focus on the 
extreme functions because they serve as strong cut-generating functions for 
general integer linear programs. (A related notion, facets, has been studied 
in parts of the literature. For the finite group problem and for special cases 
of the infinite group problem, it is known to be equivalent to that of extreme 
functions; see [9] section 2.2.4] and the forthcoming paper [27]. 1 


1.2. Extreme functions and their slopes. In this paper, we discuss how 
computer-based search can help in finding extreme functions. The next two 
subsections (1.3 and 1.4) are devoted to a short literature review on the 
success of computer-based search in these problems. 

An important statistic that has received much attention in the literature is 
the number of slopes of an extreme function. For the infinite group problem, 
we use the term k-slope function to refer to a continuous piecewise linear 
function with k different slope values, whereas for the finite group problem, 
we use the same term to refer to a discrete function whose interpolation has 
k different slope values. Figure 1 shows a 2-slope function for the finite (left) 
and infinite (right) group problems respectively. The number of slopes can 
be taken as a measure of complexity of an extreme function. Functions 


^For simplicity of notation in the n-dimensional (n-row) problem, j|5] and ESI work 
with Rf(M. n ,Z n ) instead of the aggregated formulation Rf( R n /Z n ). The aggregated for¬ 
mulation is of interest for the present paper, since for a master finite group problem where 
G = -Z, the group G/Z is indeed finite. 

"Throughout this paper, we refer to an extreme function or a family of extreme func¬ 
tions by the name of the SageMath function in the Electronic Compendium |261 i 35] that 
constructs them; these names are shown in typewriter font. The reader is invited to ex¬ 
plore these functions alongside reading this article. The Electronic Compendium is part 
of our software [23], which allows to test extremality and which has been used to make 
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Figure 1 . The 2-slope extreme function gj_2_slop^J dis¬ 
covered by Gornory and Johnson L2JJ- Left, gj_2_slope 
for the finite group problem with q = 5 and / = |, ob¬ 
tained by restrict_to_finite_group(gj_2_slope ()). It 
is a discrete function whose interpolation is the right subfig¬ 
ure. Right , gj_2_slope for the infinite group problem with 
/ = It is a continuous piecewise linear function with two 
slopes, although it has four pieces. Its restriction to is 
the left subfigure. 


with only 2 slopes are the easiest to analyze; in fact, the Gornory-Johnson 
Two Slope Theorem m states that any continuous 2-slope function that is 
minimal valid for Rf(WL/ Z) is already extreme. Gornory- Johnson in their 
2003 paper m section 6.2], after discussing the tools available for analyzing 
minimal functions for the infinite group problem, explain: 

[T]heir application to generate facets is still somewhat ad hoc. 

Also we don’t have general theorems for three or more slope 
facets analogous to the Gomory-Johnson Two Slope Theo¬ 
rem. 

Indeed, functions with 2 or 3 slopes were the focus of much of the literature, 
as we will discuss in more detail below. 


1.3. Computer-based search used in the finite group problem. Go- 

mory’s seminal paper m Appendix 5], introducing the group problem and 
corner polyhedra, listed all extreme functions up to automorphism and ho¬ 
momorphism for the finite group problems of order q = 2, 3,... ,11. Gornory 
proved that the set II^(^Z/Z) of minimal functions for R/(^ Z/Z) is a poly¬ 
tope, defined by linear inequalities that express subadditivity and certain 
equations that come from normalization; 
tails. By [T8[ Theorem 


see 


Theorem 2.2 below for de¬ 
ls]^] and |191 Theorem 2.2], the extreme functions 
for ify(^Z/Z) are the extreme points of the polytope IIy(^Z/Z). Gornory 


most of the diagrams in this paper. The captions of some figures show SageMath code 
that can be used to reproduce the diagrams. 

’hn [18] and [16] below, valid inequalities are not normalized to have the right hand 
side of 0 being 1. We state their results in our unified notation. 
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reported that the extreme points were computed, by enumerating simplex 
bases, using a computer code of Balinski and Wolfe. 

During the revival of the interest in the group problem in the 2000s, Evans 
m used her specialized implementatior|^]of the double description method 
(see, e.g., m) to enumerate all extreme points of the polytope II^(|Z/Z), 
thereby obtaining all the extreme functions for the finite group problems of 
order q < 24. By exploring the patterns of such functions, some parametric 
families of 2-slope and 3-slope extreme functions for finite group problems 
were constructed. Extreme functions from these families were generated 
by the Matlab code in m Appendix B.l] for the finite group problems of 
order q < 30. Evans reported that these extreme functions received a large 
percentage of hits in the so-called shooting experiment [161 Table 13]. 

Gornory and Johnson [19] showed that the number of extreme functions 
grows exponentially with q. Hence it is impractical to enumerate all ex¬ 
treme functions for ify(|Z,Z) when q is large. The shooting experiment 
was conducted in [22] (more results appeared in [IB] ) to identify the “im¬ 
portant” extreme functions for the finite group problems where the order is 
at most 30. This experiment was extended to finite group problems of order 
up to 90, and then to problems of order up to 200 with the so-called “par¬ 
tial shooting” variant in EO; see also [32]. Extreme functions resulting from 
the shooting experiments were expected to be important computationally in 
branch-and-cut (see, e.g., EDI Section 19.6.2] for a summary), though actual 
computational uses never seem to have materialized. They are mostly gmic 
functions (up to multiplicative homomorphism and automorphism^]) , along 
with some other 2-slope and 3-slope extreme functions. The shooting ex¬ 
periment, however, is not suitable for finding functions with many slopes, as 
those appear to be extremely rare from the viewpoint of shooting, or func¬ 
tions with specific properties for finite group problems. It is not possible to 
perform the experiment for the infinite group problem either, according to 
[30, Section 19.6.2.1], 

Araoz, Evans, Gornory and Johnson [lj demonstrated a close relation 
between the master finite group problem and the master knapsack problem. 
In particular, the convex hull P(K q -\) of solutions to the master knapsack 
problem of size g— 1 is a facet of the master finite group problem i?j(IZ/Z), 

where / = Thus, extreme functions for i?y(|Z/Z) where / = are all 
valid for the knapsack problem K q -\. Furthermore, some extreme functions 
for Rf^Jj/h) may be obtained from facets for P{K q f) through a process 
called tilting (see [U Theorem 5.2]). Examples of extreme functions derived 
by tilting a knapsack facet are listed in PQ Table B.2], 

4 lf6j Chapter 4] used a variation of the double description method that includes a 
parallel implementation for maintaining the minimal system of generators. Evans reports 
that the parallel version achieved a speedup by a factor of 12.79 using 32 processors. 

'’These two procedures are available as multiplicative_homomorphism and 
automorphism, respectively, in the accompanying SageMath program. 
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A different approach was followed by Richard, Li and Miller m, who 
proposed an approximate lifting scheme that converts certain superaddi¬ 
tive functions into potentially strong valid inequalities. The superadditive 
functions that they studied were the DPL n functions with 4n non-negative 
parameters, and the superadditive CPL n functions as a special case of DPL„ 
where 2 n parameters were fixed to 0. The parameters that define a superad¬ 
ditive DPL n function belong to a certain polyhedron PQ n (or, in the case of 
superadditive CPL„ function, to a simpler polyhedron that is a face of PQ n ). 
Several classes of well-known cutting planes can be generated by convert¬ 
ing the DPL n or CPL n functions that correspond to the extreme points of 
PQ n - However, the functions generated by the approximate lifting scheme 
are not always extreme. By the lack of any automated extremality tests for 
a parametric family of functions, the study was restricted to so small n that 
manual inspection of extremality became possible. The authors investigated 
the CPL 2 functions for the finite group problem in [29] and a special case 
of CPL 3 functions for both finite and infinite group problems in [31], all of 
which required extensive case analysis for extremality conditions by hand. 
They found the first parametric family of 4-slope extreme functions for the 
finite group problem in m- 


1.4. Computer-based search used in the infinite group problem. 

Computer-based search has also been used in the study of the infinite group 
problem. One approach focuses on continuous piecewise linear functions 
7 r with breakpoints in ^Z for some fixed q £ Z+. (Then, without loss of 

generality, one can assume that also / £ ^Z 0 Lemma 2.4].) Gomory 


and Johnson [1T9] established a connection between finite and infinite group 
problems by studying the restriction and interpolation of valid functions. 
They proved that a continuous piecewise linear function 7 r with breakpoints 
and / in ^Z is minimal if and only if its restriction to ^Z is minimal for 
the finite group problem i?j(^Z/Z). Moreover, if 7r is extreme for Rf(M./ Z), 
then the restriction 7 r|i z is extreme for i?j(^Z/Z). Therefore all search 


approaches for the finite group problem, described above in |section L3 


yield candidates for extreme functions for the infinite group problem. 

Using this connection with the finite group problem i?j(^Z/Z) amounts 
to discretizing the space of functions 7 r, by discretizing the breakpoints of 7r 

A function 7r is then uniquely determined by its 

r i —1 


in ^Z for some fixed q. 


values at ^ for i = 0 , 1 ,..., q, or by its slope values on the intervals [^-, |] 
for i = 1 ,...,q. We introduce the following notation for use throughout 
the paper. Denote the function value at | by 7 U for i £ {0,1,..., q}, where 
tq = Tj = 0. Let Si be real numbers such that qsi are the slope values on 


[■Vd for i £ { 1 ,..., q}. See Figure 2 


for an illustration. 


Chen [12] designed an enumerative algorithm to find candidate piecewise 
linear extreme functions for the infinite group problem using an additional 
discretization. In addition to fixing q, also pick a natural number v. Then 






MATTHIAS KOPPE AND YUAN ZHOU 



88888888 


Figure 2. The q x v grid discretization of the space of con¬ 
tinuous piecewise linear functions with rational data. Here 
q = 8 and v = 6. 


consider the continuous piecewise linear functions tt that have q x v grid 
discretization: the breakpoints being in {0, 1} and 

(1) TTi <E {0, i,..., 1} for i € {0, or 

(2) Si <E {-1, • • •, ^r,0, l ,..., 1} for i <E {1,... ,q}. 


In fact, these two natural ways of discretization are easily seen to be equiva¬ 


lent (Lemma B.l in Appendix B). Chen’s algorithm then enumerated every 


candidate function tt such that n is symmetric, 7r(0) = 0, n (/) = 1, 7r(l) = 0, 
and tt has the steepest positive and negative slopes at 0 and 1 respectively. 
With q = 10, v = 9, almost 500 functions were found. However, no results 
were stated regarding the extremality of these candidate functions. In fact, 
it was not until the breakthrough algorithmic results in [7] that an auto¬ 
mated test for extremality for the infinite group problem became possible. 
Despite the lack of automated test, Chen constructed the first parametric 
family of 4-slope extreme functions^ for the infinite group problem in 
Gornory-Johnson, in their 2003 paper [2L , section 6.3], had written: 

[W]e have discussed mainly large facets, families of two- 
slope facets or three-slope facets. What about much smaller 
facets? It seems extremely likely that there are much smaller 
more complex facets all around the big ones. 


However, for a brief period, it seemed that with Chen’s functions the peak 
of complexity had been reached, as no extreme functions with more than 
4 slopes were found. In 2009, Dey-Richard in m stated as an open ques¬ 
tion whether there exist extreme functions with more than 4 slopes. The 
“4-slope conjecture”, asserting that 4 slopes is the maximum for continu¬ 
ous piecewise linear extreme functions, circulated among researchers for a 


®The function is available in the Electronic Compendium as chen_4_slope 



































COMPUTER-BASED SEARCH FOR EXTREME FUNCTIONS 


9 


while. The conjecture reflected the hope that, in contrast to the “unstruc¬ 
tured and arithmetic” finite group problem, the complexity of the extreme 
functions of the infinite group problem would be under control. It is unclear 
how much support the conjecture had at any time; Dey (2016, Personal 
Communication) remembers that he did not believe in this conjecture. 

The 4-slope conjecture had to be revised quickly when new computational 
tools for testing extremality became available. Basu et al. 0 Theorem 1.5] 
(see also m Theorem 8.5]) showed that in order to test extremality of 
7 r for Rf(M./ Z), one simply needs to test extremality of its restriction to 
for Z/Z). (Later it was shown that restricting to ^ suffices ]10l 

Theorem 8.6].) This extremality test can be done by computing the rank of 
a finite-dimensional system of linear equations. 

These algorithmic results in [7] enabled Hildebrand to run a new computer- 
based search, using Matlab programs (2013, unpublished, reported in [£i 
Table 4]). Like Chen |12j . Hildebrand generated functions 7r on the q x v 
grid. However, instead of searching exhaustively, he generated the functions 
randomly. For each candidate function 7r, he checked if 7ij i z is minimal and 

extreme for Rf{^ Z/Z) (a necessary condition), and finally tested if vr|j_ z 
is extreme for Rf(^ Z/Z) (Basu et al.’s necessary and sufficient condition). 
In this way, Hildebrand discovered the first 5-slope extreme function^ thus 
refuting the 4-slope conjecture. 

Still it seemed possible that a version of this conjecture might be true; ei¬ 
ther with 4 replaced by another small number, or that the 4-slope conjecture 
holds “generically” (see [HI Open Question 2.16]). 

1.5. New search strategies; outline of the paper. In this paper, we de¬ 
velop new search strategies that aim to find extreme functions for the infinite 
group problem with many different slope values or with special properties. 
Our goals are markedly different from those of some earlier research surveyed 
above, in particular from studies using variants of the shooting experiment 
P2I221. We merely wish to settle theoretical questions regarding the struc¬ 
ture of extreme functions. We make no claims whatsoever regarding the 
computational usefulness of the functions that are found by our search. 

Our implementation is based on the software [23], which implements an 
automated extremality test, following the ideas of the proof of 0 Theorem 
1.3]. The practical implementation, described in detail in [24 , 25] , has an 
empirical running time that does not strongly depend on q, and therefore 
is suitable for functions with extremely large q , even though no theoretical 
worst-case running time bound better than polynomial in q is available. 

Like Gomory [18] and Evans m, and in contrast to Chen’s and Hilde¬ 
brand’s approaches, our strategies only discretize the breakpoints into ^Z. 

7 The functions are available in the Electronic Compendium j35] as hildebrand_5_ 
slope... 
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This is crucial to be able to reach larger values of q, as we explain in Ap¬ 


pendix B| where we discuss the complexity of search approaches based on 


q x v grid discretization. 


1.5.1. Section^ Vertex filtering search. As a first step of our study, we 
investigate how far we get by just combining state-of-the-art vertex enumer¬ 
ation software with the automated extremality test provided by the software 
[23] . We build on state-of-the-art vertex enumeration software, Parma Poly- 
hedra Library (PPL) f3| and lrslib [2j[3j. The vertex filtering search was able 
to enumerate extreme functions with q < 27, among which the first 6-slope 
extreme functions were found, breaking the previous record of 5 slopes 
due to Hildebrand (2013, unpublished; reported in [9j). From the results 
obtained by this search, we observe: 

• a diminishing fraction of vertices of Hj(^Z/Z) that correspond to 
extreme functions for Rj(R/ Z); 

• an exponential growth of time spent on vertex enumeration; 
when q increases. 


1.5.2. Section [3}’ Search using MIP. These factors suggest that one needs 
to consider other search strategies to reach larger q. We investigate search 
strategies that, for the first time, are guided directly by the subtle 
structure of minimal functions that were exposed by the proof of 
the algorithmic extremality test in [7], rather than using the extremality 
test merely as a black box. 

To this end, we review the notions of the two-dimensional polyhedral 
complex AV and of additive faces in |section 3.1 Identifying the additive 
faces of AV is a crucial step in the algorithmic extremality test 0- By 
means of Gomory-Johnson’s celebrated Interval Lemma, additive faces give 
rise to “affine-imposing” intervals (in the terminology of 0)- This reduces 
the infinite-dimensional test to a finite-dimensional one. The combinatorics 
of the additive faces of the complex AV has a central role in our new ap¬ 
proaches. 

In the remainder of Isection 3l we describe how the search based on the 
combinatorics of the additive faces can be implemented using standard MIP 
modeling techniques and running a commercial MIP solver. This is easy to 
implement and easy to tailor to a search for extreme functions with particu¬ 
lar properties. However it is limited because floating-point implementations 
are not a good match for finding functions of high arithmetic complexity, 
and because MIP solvers are generally not the best tool for performing an 
exhaustive search. 


1.5.3. Section Backtracking search. To address the shortcomings of the 
MIP approach, we have developed a specialized backtracking search algo¬ 
rithm, which we describe in section 4 Like the MIP approach, it is based 
on the combinatorics of the additive faces. It works best when combined 
with the vertex filtering search described in section 2 This combined search 
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FIGURE 3. A 28-slope extreme function kzh_28_slope_l 
found by our search code. Each color in the plotting cor¬ 
responds to a different slope value. 


algorithm looks for extreme functions by backtracking on the additive faces 
of AV in a first step and vertex enumeration in a second step. The synergy 
and balance between these two steps to obtain the best computational per¬ 
formance is discussed in section 4.7 Using the combined search algorithm, 


we discover new extreme functions with up to 7 slopes. 

1.5.4. Section^ Targeted search with patterns. In the library of functions 
found by the above algorithms, we observe some special combinatorial pat¬ 


terns on their two-dimensional polyhedral complexes AV. In section 5, we 
describe how we use these patterns to make a targeted search for functions 
with a very large number of slopes, which discovers piecewise linear extreme 
functions with up to 28 slopes. 

1.6. Summary of new results. 

Theorem 1.1. There exist continuous piecewise linear extreme functions 
with 2, 3, 4, 5, 6, 7, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, and 28 slopes. 


|Figure 3 shows one 28-slope extreme function found by our code, with 


q = 778, out of reach for any previous study. We remark that Basu et al. [6] 
have improved |Theorem 1.1| by constructing a family of extreme functions 
with an arbitrary prescribed number of slopes. 

Our computer-based search also can be tailored to find extreme functions 
with certain properties. In particular, several open questions are resolved by 
such newly discovered extreme functions. Let m > 3 be a positive integer. 
m Theorem 8.6] states that n is extreme for f?j(M/Z) if and only if the 
restriction 7r| i v is extreme for the finite group problem Rf(^~ Z/Z). Our 


i 

mq “ 


search found a function (see Figure 4) that is not extreme for Rf(W/ Z), 


but whose restriction to t^Z is extreme for Rf(^j Z/Z). This proves the 
following result, thereby answering the Open Question 8.7 in |8j. 
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Figure 4. The example kzh_2q_example_l, showing that 
an oversampling factor of m = 3 in m Theorem 8.6] is best 
possible. 


Proposition 1.2. The hypothesis m > 3 in m Theorem 8.6] is best pos¬ 
sible. The theorem does not hold for m = 2. 


The search also found piecewise linear extreme functions of Rf( R/Z) to 


answer the Open Question 2.16 in jS], which we mentioned in section 1.4 


1.7. Available software. We have made all of the discovered functions 
mentioned in this paper available as part of the Electronic Compendium 
[35]. The reader is invited to investigate the functions using our software 
[23]. The articles (24[ [25 ] describe the software in detail. The computer- 
based search code will be released shortly as part of a new version of the 
software [23]. 


2. Restriction to q grid - Vertex filtering search 

Recall that we are looking for continuous piecewise linear functions it : R — > 
M_l_ with breakpoints in |Z that are extreme for the single-row Gomory- 
Johnson infinite group problem. The construction of parametric families of 
extreme functions (a focus of many previous studies), extreme functions with 
irrational breakpoints (for example, the function 
and non-piecewise linear extreme functions such 
[5] are beyond the scope of this paper. 

2.1. Restriction to grid. Our approach is based on the discretization of 
the breakpoints of tt. More precisely, we only focus on the functions n 
with rational breakpoints in ^Z for some q E N. Suppose without loss 
of generality [TJ Lemma 2.4] that / £ ^Z. Under such hypotheses, 7r is 
uniquely determined by its values at points in ^Z. We say that 7r is the 
(continuous) interpolation of 7r|i z , while 7r|i z is the restriction of 7r to the 


Q Q 


grid iZ. 

Figure 1 

in 

section 1 

illustrates the interpolation and restriction of 

a gj_2_slope function with g 

= 5. 


Gornory and Johnson proved the following relations between 7r and 7r|i z : 


bhk_irrational in 0), 
as bccz_counterexample 
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Theorem 2.1 (p9|; see also m Theorem 8.3]). Let n be a continuous 
piecewise linear function with breakpoints in f Z for some q E Z + and let 
f E |Z. Then the following hold: 

(1) 7 r is minimal for Rf(M./Ti) if and only if 7 r|i z is minimal for i?/(-Z/Z). 

(2) If it is extreme for i?r(M/Z), i/ien7r|i z is extreme for i?/-(-Z/Z). 

Hence the interpolations of those 7 r| i z that are extreme for Z/Z) are 

the only possible candidates of continuous piecewise linear extreme functions 
with breakpoints in ^Z for ify(M/Z). Extreme functions are clearly minimal 
functions; the latter have a characterization given by a theorem by Gomory 
and Johnson. In the case of the finite group problem Rf(^ Z/Z), it can be 
stated as follows. 


Theorem 2.2 (|19|; see also |9, Theorem 2.6]). Let it and f be as above. 
7r|i z is minimal for Rj(J^L/Tl) if and only if 

(1) 7T 0 = 0, 

(2) 7r| i z is subadditive: n^ x+y ) mod q < n x + n y for x, y E Z, 

(3) 7r| i z is symmetric: tt x + 7 T q f- X = 1 for i£Z, 

9 

where 7r* = 7r(|) for i E Z. 

Since 7r: M/Z —>• M is periodic modulo 1, a minimal function 7 r|i z for the 

9 

finite group problem is specified by its values (7 To,7Ti, ..., 7t ? _l) on the grid 
points ^Zn [0,1). The following statement immediately follows from the 
observation that the above conditions are all linear constraints. 


Proposition 2.3 ( |19l Theorem 2.2]). The set n^QZ/Z) of minimal func¬ 
tions for Rf(^j Z/Z) is a convex polytope. Furthermore, extreme functions 
for Rf(^ Z/Z) are the extreme points (i.e., vertices) of this polytope. 


By |Theorem 2. 1| and [Proposition 2.3 , all continuous piecewise linear ex¬ 
treme functions for Rj(R/ Z) with breakpoints in ^Z can be found by inter¬ 


polating the vertices of the polytope n^(^Z/Z). However, in general, ex- 
tremality of 7 r| i z for (f Z/Z) does not imply extremality of 7 r for Rf(R/ Z). 

a " 


This makes further hltering necessary, which we explain below. 


2.2. Preprocessing. We begin with a remark on the minimal H-represen- 
tation of n^(^Z/Z). The H-representation of the polytope nj(-Z/Z) de- 

l. 


fined by 


Theorem 2.2 


has asymptotically fq 2 constraints, many of which 


are redundant, indeed, m Corollary 2.7] gives a minimal representation of 
n f (±Z/Z) that only has asymptotically |g 2 constraints, mainly by replacing 


the subadditivity constraints ([2]) of Theorem 2.2 


Hi + 7 Tj > 7r (i+i ) mod g for o <i<j<q 
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with the triple system: 


7Tj + TTj + TTk > 1 for 0 < i < j < k < q, i + j + k = qf (mod q). 

A minimal H-representation is of interest for vertex enumeration, because 
having many redundant inequalities may greatly slow down the vertex enu¬ 
meration process. Although a minimal H-representation is known for the 
polytope n^(^Z/Z), the search strategies described in later sections of the 
present paper also need to deal with other polytopes whose minimal H- 
representations are not known. For this purpose, computational preprocess¬ 
ing is used to remove the redundant inequalities from the H-representation 
of a polytope before enumerating its vertices. Namely, we apply the pre¬ 
processing program redund provided by lrslib (version 5.(Q, which removes 
redundant inequalities using Linear Programming. The third columns of Ta¬ 
bles [l] and Table 2| show the number of inequalities in the H-representation 
of the polytope ny(|Z/Z) for / = ^ before and after preprocessing, respec¬ 
tively. The number after preprocessing is roughly i q 2 , which is consistent 
with m Corollary 2.7], 


2.3. Performance of various vertex enumeration codes. Various soft¬ 
ware packages are available for computing the vertices of a polytope given 
by an H-representation. We considered the following popular packages. 

• cddlib (version 094^]), a well-known implementation of the double 

description method M- „ 

• Parma Polyhedra Library (version l.lp^|). PPL is a C++ library 
for the manipulation and computation of rational convex polyhedra 
gj. Like cddlib, polyhedral computations in PPL are based on the 
double description method. 

An extensive computational study g], section 4] showed that the 
double description method implementation in PPL has a better per¬ 
formance (on the vertex/facet enumeration problem) compared with 
cddlib and with other polyhedra libraries that are popular in PPL’s 
primary application domain, New Polka and PolyLib. 

• Porta (version f ,4.1p^|, based on the Fourier-Motzkin elimination 
method. 

• lrslib (version 5.0), a C implementation of the lexicographic reverse 
search algorithm for vertex enumeration and convex hull problems 
BIG]. This algorithm uses little memory space during the com¬ 
putation; vertices are generated as a stream and are not stored in 
memory, which makes it suitable for vertex enumeration problems 
for polytopes with a large number of vertices. 

® Available from http: //cgm. cs .mcgill. ca/~avis/C/lrs .html. 

® Available from http://www.inf . ethz . ch/personal/fukudak/cdd_home/cdd.html 

^Available from http://bugseng.com/products/ppl/ 

11 Available from http: //www. iwr.uni-heidelberg.de/groups/comopt/software/ 
PORTA/ 
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The computational study [H, Table 2] showed that lrslib outper¬ 
forms PPL for large problems, whereas PPL outperforms lrslib for 
easy problems jl; Table 1]. 

Based on the study [3], we expected PPL to be the best choice for low 
dimensions and lrslib to be the best choice for high dimensions. We decided 
to verify this using a computational study for our polytopes L1 j(3z/Z) given 
Proposition 2.3 for various values of q and / = -, recording the running 


m 


times for vertex enumeration. 

In addition to the four packages listed above, we also included the follow¬ 
ing package in our experiments, which was developed very recently and had 
not been investigated in a major computational study. 

• PANDA (version 2015-02-2^^) based on the parallel adjacency de¬ 
composition algorithm 


At the time of the revision of the present paper, we also became aware of 
the high performance vertex enumeration code in the following package. 

• Normaliz (version 3.1.1p*|), a software tool for the computation of 
Hilbert bases and enumerative data of rational cones and affine 
monoids. Its vertex enumeration code is based on the Fourier- 
Motzkin elimination method and the pyramid decomposition de¬ 
scribed in mi- 


We used PPL via its Cython-based library interface within the SageMath 
computer algebra system. The other systems were tested using their 
command-line executables because no library interfaces were available for 
them in SageMath j 71 ] Hence in very low dimension, there is a slight bias in 
favor of PPL due to the overhead in using the other systems via executables 
and file passing. All systems were tested using exact computation modes 
and using a single thread. 

Tables |T] and [2] report for each test the size of the polytope and the run¬ 
ning times measured in CPU secondjj^J without and with preprocessing, 
respectively. The preprocessing in |Table 2| consists of removing redundant 
inequalities from the H-representation using the command redund provided 
by lrslib. We also measured the computational overhead of interfacing to 
redund in Python, which exceeds the actual redund running times. Compar¬ 
ing Table 1 and Table 2 shows that for dimension at most 6, it is best to just 
run PPL on the original inequalities without preprocessing. For dimensions 


12 Available from 

1 ^ Available from 

14 A dataset with the input hies in the formats of the systems in our study is 
available at https : //www.math .ucdavis . edu/~mkoeppe/art/inf inite-group/dataset_ 
gomory_polyhedra_cyclic_group.zip 

10 The tests have been performed on a virtual machine running under the QEMU hy¬ 
pervisor, which reports to have access to 12 processors running at 2.0 GHz. However, 
due to the virtualization, the measured running times have a large variance between runs, 
though all algorithms are deterministic. 


http://comopt.ifi.uni-heidelberg.de/software/PANDA/ 
https://www.normaliz.uni-osnabrueck.de/ 
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at least 7, Normaliz is the clear winner. The total time for preprocessing 
pays off when the dimension of the polytope is at least 10 for Normaliz, or 
is at least 8 for the other systems. For dimension at least 13, Normaliz is 
faster than the second-best code, lrslib, by more than an order of magnitude. 
Normaliz is also the least dependent on preprocessing!^] 

For the computational experiments in the remainder of this paper, we 
used a combination of PPL (for low dimensions) and lrslib (for high dimen¬ 
sions), using preprocessing when the dimension is at least 8. We did not use 
Normaliz, as we only became aware of its high performance code for vertex 
enumeration at the time of revision of the present paper. 


2.4. Filtering. As mentioned earlier, extremality of vr|i z for i?*(-Z/Z) 
does not always imply extremality of n for ify(M/Z). We call the interpola¬ 
tion 7r of a vertex ir\ i z of the polytope n^(^Z/Z) a vertex-function. Once the 

vertices of Flj(^Z/Z) are enumerated, we can use the automated extremality 
test implemented in the software [231425] to filter out those 7 r|i z whose in¬ 
terpolations are not extreme for itb(M/Z). We will discuss in 


section 3.1 


the 

two-dimensional polyhedral complex AV and the notion of covered intervals 
that make this automated extremality test possible. 

Remark 2.4. Given that 7r| i z is a vertex of IIj(4z/Z), the extremality test 

of 7r for i?j(M/Z) reduces to testing whether all intervals are covered, accord¬ 
ing to |Theorem 3.2 This can be determined using the function generate_ 
uncovered_intervals of 


We now summarize the above ideas in the following algorithm, which is 
referred to as “vertex filtering mode” in our code. The implementation uses 


Parma Polyhedra Library and lrslib as described in section A.l 


Algorithm 1: vertex filtering mode 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


Consider the restriction of 7r to the grid 


Define 7To, 7 Ti, ..., ir q as variables, where 7u = 7 r(^). 
Construct the polytop e IIj(4z/Z) o f minimal functions for 
i?y(4z/Z), defined by 


Enumerate the vertices 7r|i 

9 

For every vertex 7 r|i y , do: 


Theorem 2.2 

of this polytope. 


(a) Interpolate to get 7r, a minimal valid function for i?y(M/Z). 

(b) If the intervals [|, for i = 0,1,..., q — 1 are all covered, 
then 7r is extreme for i?y(M/Z). Output function 7r. 


Bruns (Personal Communication, 2016) explains that Normaliz performs initial 
transformations that remove a large part of the redundancy present in these problems. 
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Figure 5. Vertex enumeration time (not including checking 
extremality of vertex-functions), mean extremality checking 
time for a vertex-function, and percentage of extreme func¬ 
tions 


2.5. Performance of the vertex filtering search. Our search code is 
implemented in SageMath [jjPQ, an open-source mathematics software system 
that uses Python and Cython as its primary programming languages and 
interfaces with various existing packages. 

Our vertex filtering search code uses the strategies described in 
to decide whether preprocessing is needed and which software 
vertex enumeration. We test its performance for q = 10,11,..., 27 and 


section 2.3 
to use for 


/ = | for x = 1,2,..., UJ. 

Observe that as q increases, the dimension and the number of vertices 
of the polytope increase. In particular, it results in an exponential growth 
of running time for vertex enumeration (cf. Figure 5-left). In addition to 
vertex enumeration, the vertex filtering search has to run extremality tests 
for the vertex-functions once they are found, which consumes non-negligible 
extra time (cf. Figure 5-middle). Furthermore, Figure 5-right illustrates 


a decrease in the percentage of extreme functions to vertex-functions. It 
suggests that when q is large, vertex filtering search does enumeration in 
high dimension and throws away many non-extreme functions. Therefore, it 
is not surprising that the vertex filtering search is only suitable for small q 
(q < 27). 


2.6. Results and discussion. Despite its limitations, vertex filtering search 
finds up to 6-slope extreme functions with q < 27, breaking the previous 
record of 5 slope?) 1 7 | due to Hildebrand (2013, unpublished). 

The extreme functions obtained from the vertex filtering search with 
breakpoints in ^Z/Z (q < 27) for the infinite and finite group problems 
are plotted in |Figurc 6| using blue and red dots. They are placed from left 
to right according to the value q and the number of slopes. The log-scale 
y-axis refers to the arithmetic complexity of the extreme functions, which 
we define as follows. 


^These functions are available in [33] as hildebrand_5_slope. . . 
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Definition 2.5. The arithmetic a function ir: ^Z/Z — > [0,1] 

(or of a continuous piecewise linear function n: M/Z —> [0,1] with break¬ 
points in ^Z) is defined as the least common denominator of the values 
TTi = 7 r(^) for i E {0,..., q}. 

For example, it is easy to see that the gmic function with / E ^Z has an 
arithmetic complexity of 0(q 2 ). 


We observe from Figure 6 an empirical exponential growth of the arith¬ 


metic complexity as q increases. This suggests that huge values of v would 
be needed in a search based on the q x v grid discretization, like Chen’s and 


Hildebrand’s described in section 1.4), in order to make them exhaustive. 


This renders this approach unsuitable for large q. 

In the following, we give an upper bound on the arithmetic complexity of 
an extreme function in terms of q. We will refer to this boun d in the following 
sections. Let 7r be an extreme function for Rf(^ Z/Z). Proposition 2.3 

states that (7 To,7Ti, ... ,7 r g _i) is a vertex of the polytope defined 

in | Theorem 2.2| By introducing slack variables, the constraint system of 
n/(iZ/Z) can be written in the standard form using matrix notation as 
Ax = b, x > 0, where A and b have all integer entries. Then by Cramer’s 
rule, the denominators of {vr^l come from the inverse of simplex 
basis matrices. The next lemma investigates the determinants of simplex 
basis matrices of A. The proof of this lemma appears in Appendix B| 

Lemma 2.6. Let q E Z + and f E ^Z, 0 < / < 1. Let Ax = b, x > 0 be 
the constraint system of \Theorem 2.2 \ written in the standard form. Let B 
be a basis matrix of A. Then |detL>| < IO 9 / 4 . 

Given q E Z + , let d veT and d ex t denote the maximum arithmetic com¬ 
plexities of any extreme function 7r for Rf(^j Z/Z) and of any extreme func¬ 
tion 7r for Rf(M./ Z) with breakpoints in ^Z, respectively. By 


Theorem 2.1 

and Proposition 2.3] it is clear that the number d ex t is well-defined and 
dext < d ver . Since the entries of the right-hand side b of Ax = b are inte¬ 
gers, we have that 

dext < d V er < dbas := max{ |det B\ : B is a basis matrix of A } < IO 9 / 4 . 
We state this upper bound on d vei and on d ext as a corollary. 

Corollary 2.7. Let it be an extreme function for Rf(M./ Z) with breakpoints 
in ^Z (or an extreme function for Rf(^h/h)). Then the arithmetic com¬ 
plexity of ir (i.e., the least common denominator of (ttq, tt\, ..., n q )) is at 
most IO 9 / 4 . 


We do not have a matching lower bound for d ex t; however, Lemma B.2 


m 


Appendix B| offers an exponential lower bound on dbas when an arithmetic 
condition on q and / is satisfied. 


^This function is available as arithmetic_complexity 

























Table 1. Efficiency of various vertex enumeration codes without preprocessing 


q 

dimension 

inequalities 

vertices 



Running time (s) 



PPL 

Porta 

cddlib 

lrslib 

Panda 

Normaliz 

5 

1 

21 

2 

0.001 

0.018 

0.009 

0.008 

0.026 

0.003 

7 

2 

36 

4 

0.001 

0.012 

0.011 

0.005 

0.026 

0.004 

9 

3 

55 

7 

0.002 

0.016 

0.018 

0.004 

0.065 

0.005 

11 

4 

78 

18 

0.003 

0.016 

0.031 

0.009 

23 

0.007 

13 

5 

105 

40 

0.007 

0.018 

0.11 

0.021 

4604 

0.011 

15 

6 

136 

68 

0.017 

0.037 

0.21 

0.14 


0.017 

17 

7 

171 

251 

0.14 

0.20 

1.2 

0.71 


0.047 

19 

8 

210 

726 

0.91 

1.6 

5.0 

2.3 


0.16 

21 

9 

253 

1661 

6.6 

13 

24 

13 


0.67 

2m 

10 

300 

7188 

166 

558 

785 

74 


4.9 

25 

11 

351 

23214 

1854 

10048 

12129 

471 


21 

26 

12 

378 

54010 




2167 


62 

27 

12 

406 

68216 






89 

28 

13 

435 

195229 






326 

29 

13 

465 

317145 






644 

30 

14 

496 

576696 






1693 

31 

14 

528 

1216944 






3411 


“By isomorphism, this vertex enumeration problem (q = 23, / = is the same as the problem with q = 23 and f = The latter was tested by 
L. Evans ED Table 6] using her parallel C implementation of the double description method, reporting a running time of 9.58 hours (ca. 34500 s) 
on one processor and 0.75 hours on 32 processors, each a 550MHz Pentium III Xeon, on the Jedi cluster of the Interactive High Performance 
Computing Cluster at Georgia Tech. 
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Table 2. Efficiency of various vertex enumeration codes with preprocessing 


to 

o 


Running time (s) 


q 

dim 

ineqs 

vertices 

PPL 

Porta 

cddlib 

lrslib 

Panda 

Normaliz 

redund 

overhead 

5 

1 

7 

2 

0.003 

0.009 

0.006 

0.010 

0.019 

0.003 

0.006 

0.040 

7 

2 

10 

4 

0.001 

0.010 

0.009 

0.007 

0.015 

0.003 

0.006 

0.029 

9 

3 

14 

7 

0.001 

0.008 

0.009 

0.008 

0.021 

0.004 

0.009 

0.010 

11 

4 

20 

18 

0.002 

0.008 

0.015 

0.010 

0.017 

0.006 

0.012 

0.049 

13 

5 

27 

40 

0.003 

0.007 

0.021 

0.012 

0.039 

0.009 

0.022 

0.050 

15 

6 

35 

68 

0.004 

0.012 

0.032 

0.025 

0.040 

0.012 

0.041 

0.14 

17 

7 

45 

251 

0.016 

0.030 

0.22 

0.10 

0.16 

0.029 

0.041 

0.21 

19 

8 

56 

726 

0.061 

0.087 

0.34 

0.48 

0.44 

0.085 

0.16 

0.40 

21 

9 

68 

1661 

0.25 

0.25 

1.1 

2.5 

3.1 

0.22 

0.25 

0.72 

23 

10 

82 

7188 

4.0 

4.1 

8.0 

15 

9.0 

1.3 

0.46 

1.1 

25 

11 

97 

23214 

69 

43 

31 

94 

15 h 

6.0 

0.75 

1.8 

26 

12 

115 

54010 

511 

350 

692 

594 


21 

0.95 

2.6 

27 

12 

113 

68216 

433 

493 

672 

543 


23 

1.0 

3.0 

28 

13 

133 

195229 

8399 

5796 

9550 

3617 


102 

1.6 

4.0 

29 

13 

131 

317145 

18361 

11341 


3366 


158 

1.9 

4.9 

30 

14 

152 

576696 

> 1 d 

66747 


22743 


488 

2.5 

6.1 

31 

14 

150 

1216944 

> 3 d 

> 2 d 


20407 


734 

2.8 

7.8 
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TABLE 3. The arithmetic complexity of the search based on q x v grid discretization. d ext and d vel are the empirical 


values of v for infinite and finite group problems, and |detB| is the estimated value from Lemma B.2 


q 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

^ext 

21 

30 

35 

48 

51 

64 

63 

120 

91 

168 

165 

208 

255 

348 

289 

504 

459 

800 

dver 

21 

30 

35 

48 

51 

70 

65 

138 

95 

210 

165 

250 

315 

570 

425 

768 

651 

1120 

det-B 


32 


64 


128 


256 


512 


1024 


2048 


4096 


8192 



FIGURE 6. Arithmetic complexity and number of slopes depending on q. Extreme functions tt with breakpoints in 
^Z/Z for Rf(M./ Z) and for R/{^ Z/Z) are plotted in dark blue and bright red , respectively. The i-axis refers to the 
value q. Within the same value q, extreme functions tt are placed in ascending order by their number of slopes, from 
left (2 slopes) to right. The log-scale y-axis refers to the arithmetic complexity of the extreme functions tt, i.e., the 
least common denominator of {7To, 7 Ti, ..., 7r g _i}, showing the complexity of an exhaustive search based on q x v grid 
discretization. 
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3. MIP APPROACH 


In this section we present a new approach for computer-based search 
for extreme functions. It uses standard mixed integer linear programming 
(MIP) modeling techniques to obtain a MIP that mimics the algorithmic 
extremality test due to Basu et al. [7]. 


3.1. The two-dimensional polyhedral complex A7 7 . We first review 
the notion of a two-dimensional polyhedral complex, which serves as a tool 
for studying additivity relations and covered (affine imposing) intervals of 
piecewise linear functions. We follow [5) Section 3], but define the notions in 
our case where the function n is continuous, piecewise linear and has all its 
breakpoints in ^Z. This matches the setting of [7]. Since a minimal function 
is periodic modulo 1, we can restrict the study to the domain [0,1] only. We 
define the evenly spaced one-dimensional polyhedral complex Vi z to be the 

9 

collection of singletons and elementary closed intervals on the grid -Z by 



For any /, J, K 6 Vi z , let 

9 

F(I, J, I <) := { (x, y) G I X J : x © y G K} C [0,1] X [0,1], 


where x © y = (x + y) mod 1. Then the set 


AVi v := 


{f(I,J,K) :I,J,K€Vi z ] 


is a two-dimensional polyhedral complex. It is the collection of the ele¬ 
mentary upper or lower triangles on the grid -Z x -Z with the vertices 
(zero-dimensional faces) and edges (one-dimensional faces) that arise as in¬ 
tersections of these triangles (two-dimensional faces). See Figure 7 for an 
illustration. 

Define the subadditivity slack An: [0,1] x [0,1] —> R of n by 


A-7 t(x, y) := n(x) + n{y) — n{x © y) 


for x,y G [0,1]. Note that An is non-negative if n is minimal, since min¬ 
imality implies subadditivity. A face F of the two-dimensional polyhedral 

complex AVi z is said to be additive if An = 0 on F. Together the addi- 

q 

tive faces form the additivity domain of the function n. Since n is linear 
on the intervals of Vi z , the function An is linear on the faces of the two- 

q 

dimensional complex AVi z . Hence, the condition above is equivalent to 


Figure 7 


we m- 


An = 0 on the set of vertices of F. 

In the diagrams of the polyhedral complex AV\ Z such as 

9 _ 

dicate additive faces by various colors. Isolated additive points and additive 
edges are always drawn in green; two-dimensional additive faces (triangles) 
are painted in various colors, the significance of which we shall explain be¬ 
low. If a triangle is white, this means that strict subadditivity holds in the 
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Figure 7. Diagram of a minimal valid function ( graphs 
on the top and the left ) on the grid j^Z and the corre¬ 
sponding painting on the two-dimensional polyhedral com¬ 
plex AVj_ z ( gray solid lines), as plotted by the com¬ 
mand plot_2d_diagram(7r, colorful=True), where n = 

not_extreme_l () . Faces of AVj^ v on which A it = 0, i.e., 

io' 2 ' 

additivity holds, are shaded in colors that correspond to the 
4 connected components of this function. The heavy diagonal 
green lines x + y = f and x + y = 1 + / correspond to the 
symmetry condition. At the borders, the projections Pi{F) 
of two-dimensional additive faces are shown as gray shadows : 
Pi(F) at the top border, P 2 (F) at the left border, ps(F) at 
the bottom and the right borders. 


interior of this face. We will often refer to the diagram of the additive faces 

as the painting of vr on the complex AVi z . 

q 

An additive face implies, among other things, the important covering 
(affine imposing in the terminology of [7]) property that we outline here. 
We refer the reader to |9j, Section 4] for details on the Interval Lemma and 
its generalizations. 

Define the projections p\ ,p 2 ,P 3 '■ [0,1] x [0,1] -* [0,1] by 
Pi{x,y) = x, p 2 (x,y) = y, p 3 (x,y) = x®y. 

Let F be a two-dimensional additive face of AVi% (i.e., F is an elementary 
upper or lower triangle in the two-dimensional polyhedral complex such that 
Avr = 0 on F). By the convex additivity domain lemma for continuous func¬ 
tions B Corollary 4.9] (a consequence of the celebrated Gomory-Johnson 
Interval Lemma), ir is affine imposing with the same slope on the projec¬ 
tion intervals p\ (F) , p 2 (F) and p 3 (F). We say that these three intervals are 
(directly) covered and connected to each other. 
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Let F be a one-dimensional additive face of AV 1 z (i.e., F is an elemen- 
tary horizontal, vertical or diagonal edge in the two-dimensional polyhedral 
complex such that Avr = 0 on F). Then two of the projections Pi(F),p 2 (F) 
and ps(F) are one-dimensional. These two intervals are said to be connected 
by an edge. An interval that is connected to a covered interval is also said 
to be (indirectly) covered ; see [7J Lemma 4.5]. 

The covered intervals of n are computed in two steps. Start with di¬ 
rectly covered intervals as pi(F),p 2 (F) and ps(F) of two-dimensional addi¬ 
tive faces F. Then continue transferring indirectly covered properties using 
one-dimensional additive faces until no new covered intervals are found. 
(This saturation process clearly ends after a finite number of steps.) 

The set of covered intervals is partitioned into connected component^ \ 
In the diagrams of the polyhedral complex AVi z , colors are used to indicate 

9 _ 

membership in a connected component. The function in Figure 7 for ex¬ 
ample, has 4 connected components, though it only has 3 slopes. Within a 
connected component, the function it has the same slopes. Thus the number 
of connected components gives an upper bound on the number of slopes of 
the function 


Remark 3.1. Though the number of different slopes of a function has at¬ 
tracted the attention in the past, it appears that the number of connected 
components is a more fundamental notion. 


A painting is called a covering painting if all intervals [^, ^i] (0 < x < 
q — 1) are covered. By Theorem 3.2 below, every extreme function n has a 
covering painting. This property will be used as an important ingredient in 
the MIP approach described in this section and also the backtracking search 
approach to be discussed in |section 4 


Theorem 3.2 (rephrased from results in [7]). Let it be a continuous piece- 
wise linear function with breakpoints in ^Z for some q £ Z+ and let f € |Z. 
Then ir is extreme for Rt(M./ Z) if and only if tt\ i z is extreme for itb(-Z/Z) 

q 1 

and all intervals [|, ^±i] for x = 0,1 ,..., q — 1 are covered. 


Instead of giving a proof, we point the reader to the results from [7] that 
are rephrased as Theorem 3.2 The “if” direction follows directly from 0 


Corollary 3.4]. If vr|i z is not extreme for Rf(^ Z/Z), then 7r is not extreme 
for R f (R/ Z) by 


Theorem 2.1 


If the intervals [|, ^±1] for x = 0,1,..., q — 1 
are not all covered, then [7j Lemma 4.8] implies the nonextremality of n by 
equivariant perturbation. This shows the “only if” direction by contraposi¬ 
tion. 


^The connected components are understood in a graph-theoretic sense. We refer the 
reader to [7j for details on the graph of intervals that is used. 

■^Available as immber_of.components and number.of.slopes in [231 . 
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Px,y +1 Px+l,y+l 



Px,y h X) y Px+l,y 


Figure 8. Binary variables for color: additive (colored dark 
green in the diagram) = 0; strictly subadditive (colored white 
in the diagram) = 1 


3.2. Additivity variables and prescribed partial paintings. In the 

MIP approach we use binary variables to control the additivity (coloredness) 
of a face, i.e., a triangle, an edge, or a vertex of the two-dimensional complex 
APi z . Let x, y be integers between 0 and q — 1. We use variables k x . y for the 

lower triangle whose lower left corner is vertex (|,|), v X:V for its vertical 
edge, h x>y for its horizontal edge, d X)V for its diagonal edge, u XjV for the 
upper triangle with the same diagonal edge, and p X)V for the vertex ( x,y). 
The value 0 for these variables represents additivity (colored face) and the 
value 1 represents strict subadditivity in the relative interior of the face 
(white face); see Figure 8. These variables are subject to the invariance of 


the subadditivity slack under exchanging x and y , hence we have 


Px,y — Py,xi @-x,y — y,xi u x,y — ^y,xj hx, y — Vy t X and d Xt y — dy^ x , (5) 

for x, y G {0,1,..., q — 1}. Therefore, it suffices to consider only the upper 
left triangular part of the complex APi z where x < y. These binary vari- 


9 

ables are also subject to inclusion constraints: for any x,y € {0,1,... ,q— 1}, 
maxjpj y, p x y j r \, p x ^.\ y } ' ^x,y A Px,y T Px,y+1 T Px-\-l,yi (6a) 

max-^Pjjy.)-!, Px+1 ,y+l ) Px-\-l,y } A ^x,y A Px,y+1 T Px+l,y+1 T Px+l,yi (6b) 

inax{/) :ry . Px.y ■ 1 } A Vx,y A Px,y T Px,y+li (6c) 

max'l/Vy. p x ■ i.;y} A h :l: y < Px,y T Px+i,yi (6d) 

y . i. f) x . l.y } A d x ^ y A Px,y+1 T Px+l,yi (6e) 

where p x , q = p 0tQ , p (hy = p qfi and p q , q = p 0 , o- 


Remark 3.3. The subproblems obtained from branching on the values of 
these binary variables have an obvious interpretation in terms of paintings 
on the two-dimensional complex AVi%, which we shall refer to as prescribed 

9 

partial paintings : If an additivity variable is not fixed yet in a branching 
node, the corresponding face is shown in light grey. If it has been fixed to 0, 
the face is painted dark green. If it has been fixed to 1, the face is painted 
white. 

We shall say that a function satisfies a prescribed partial painting when 
it would be a feasible solution to the corresponding node subproblem, i.e., if 
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it satisfies all additivity conditions and strict subadditivity conditions corre¬ 
sponding to faces that have been painted dark green and white, respectively. 


Example 3.4. |Figure 10 shows a tree of prescribed partial paintings of 
APi z with q = 4 and f = \, obtained from branching on the variables £ Xjy 


and u x , y . In (the upper left triangular part of) the partial painting at its 
root node (a), the values of the following binary variables are set to 0. 


Po,o, Po,i, 


P0,2 > P0,3, P 2 , 31 ^ 0 , 0 ) ^ 0 , 1 ) '^ 0 , 2 , 110 , 3 , ^ 0 , 0 , ^ 1 , 3 , 


d-2,2- 


3.3. Function value variables. The values of candidate functions are 
modeled by continuous variables ttq, tti, ..., 7r g _i E [0,1]. They are subject 
to the constraints in IThcorem 2.21 


7Tq = 0 (7a) 

e Px,y £: Te + TTy — TT( X _|_y) mo( j g < 2p x ,y (”7b) 

Te T mod q (^) 


where e is a small positive number used to enforce the strict subadditivity 
of vertices (x,y) with p x>y = 1. See section 3.6 for a further discussion on 
the value of e. 


3.4. Slope value variables and assignment. In a search for functions 
with a prescribed number k of different slopes, we introduce k continuous 
variables si, S 2 , ■ ■ ■, Sfc E [—q, q] for the different slope values of ir. We can 
enforce si > S 2 > ■ ■ ■ > Sk by another artificial lower bound e' (see 


sec¬ 


tion 3.6): 


— 8 j+\ > e 7 , for 1 < j < k — 1. 


(8a) 


Then binary variables 6 x j (0 < x < q — 1,1 < j < k) are used to assign 
intervals to slope values: 


k 

6 x .j = 1, for 0 < x < q — 1, (8b) 

j= 1 


Sj = q(7r x+ i — tt x ) if and only if 5 x j = 1, for 0 < x < g — 1 and 1 < j < k. 
The last condition can be written as the linear inequality 

| Sj + q(n x - 7rm+i)| < 2q(l - S x j). (8c) 

We add the linear constraints 

9-1 

yy 5 x j >1, for 1 < j < k (8d) 

x=0 

to the MIP formulation, so as to ensure that every slope value Sj is used by 
at least one interval. 
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3.5. Variables for directly and indirectly covered intervals. We use 

binary variables c z ,o (0 < z < q — 1) to control whether the interval [|, 
is directly covered or not: 0 for covered and 1 for uncovered. They are 
subject to combinatorial conditions of being directly covered by dark green 
triangles presented in section 3.1 For 0 < z < q — 1, we have that c Z) o = 0 


if and only if at least one of the following binary variables has value 0. 


4 , 2 / (0 < y < q - 1), u z<y (0 < y < q - 1), (9a) 

4,z (0 < x < q — 1), u x>z (0 < x < q - 1), (9b) 

4,i/ (0 < x, y < q — 1 such that x + y = z (mod q)), (9c) 

u x , y (0 < x, y < q — 1 such that x + y + \ = z (mod q)) (9d) 


We assume that the saturation process of transferring indirectly covered 
properties described in section 3.1 ends in a finite number maxstep of steps. 
(Though no theoretical bound better than maxstep < q is known, in practice 
a small value of maxstep such as 2 is sufficient.) For 1 < i < maxstep, we 
define binary variables c z ^ (0 < z < q — 1) to model whether the interval 
[|, is covered in the first i steps of the saturation process. If the interval 

[|, was already covered in the step i — 1, then it remains covered in the 

step i. Thus, we have the constraint 


i — 0 ^ c Zt i — 0. (10a) 

If the interval [|, is connected by a green edge to an interval [|, 
that was already covered in the step i — 1, then the interval [|, becomes 
covered in the step i. 


lliill{d x z , h x ^ z _ x j mod qi V(x— z) mod q,z} 0 and C Xt i— 1 0 =>■ 0. 

(10b) 

Otherwise c Zj j = 1. 

Note that these constraints can all be expressed using linear equations or 
inequalities over binary variables. 


3.6. Trade-off between strictness and basicness. Assume that the vari¬ 
ables and the constraints defined in lsection 3.41 were not introduced to our 
MIP, and that the strict subadditivity constraints were not enforced (i.e., 
set e = 0 in ([7b]) ). 


Remark 3.5. Painting faces dark green in APi z (i.e., setting binary vari- 

ables p X)V to 0) amounts to restricting the corresponding subadditivity in¬ 
equalities ([2]) of Theorem 2.2 to equations in the constraint system of the 
polytope IIj(^Z/Z) of minimal functions. Thus, the set of restricted func¬ 
tions 7ij 1 z that satisfy the new constraint system is a smaller polytope, 


which is a face of the polytope IIj(^Z/Z). A vertex of the smaller polytope 
is also a vertex of the polytope IIj(^Z/Z). 
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After fixing 0/1 variables, a basic feasible solution ir of the system with 
c 2 ,maxstep = 0 for 0 < 2 < q — 1 is extreme for Rf{R/Z), according to 
Theorem 3.2 Unfortunately, in this case we cannot expect the solutions 


returned by the MIP solver to always have k different slope values; indeed, 
they often degenerate to 2-slope or 3-slope functions. The same phenomenon 
was observed in the shooting experiments in the literature |16l JJ/221, where 
the extreme functions that received a large percentage of hits are 2-slope or 
3-slope functions. 

Because of this, we use the slope value variables and constraints from 


sec¬ 


tion 3.4 in the MIP, always enforcing the strict inequality in slope variables 


by making a practical choice of e' > 0 for (8a). In this way we can set the 


number of slopes k of the resulting functions explicitly. We also use e > 0 


in (7b), which allows for a significant speedup. 


Remark 3.6. Corollary 2.7 shows that we do not lose any extreme functions 
by setting 0 < e, e' < 10~ 9 / 4 . Our code instead uses the heuristic choice 
e = e r = 1/4 (or 1/12), which allows for stronger pruning based on linear 
programming at the expense of losing some functions. 

However, now a basic feasible solution of the system after fixing 0/1 vari¬ 
ables is no longer guaranteed to be an extreme function, since the smaller 
polytope with the inequalities n x + ir y — 717 x + y ) mod q > e is not a face °f 
the polytope Hj(fZ/Z). In this case, a further extremality test needs to be 
applied to the returned candidate function ir. Since all intervals are covered 
(i.e., c^jnaxstep = 0 for 0 < 2 : < q — 1), it suffices by Theorem 3.2 to test 
whether vr|i z is a vertex of Hj(fZ/Z). 


3.7. Objective function. A tailored objective function can be used to steer 
the optimum away from equality of slopes and from the lower bounds of the 
subadditivity slacks An Xty , ensuring that the basic optimal solution returned 
by the MIP solver will correspond to an extreme function. However, there 
is no a priori best choice of such an objective function that will guarantee 
success. In our computations, we always maximize the difference of slopes 
si — Sfc. Other objective functions, including the following, are plausible: 


some 


fc-l 

maximize a weighted difference of slopes E A j(sj - Sj+i), for 

3= 1 

suitable weights \j; 

maximize a weighted sum of subadditivity slacks E 

0<x<y<q—l 

maximize a weighted sum of subadditivity points E ^x,yPx,yi 

0<x<y<q —1 

minimize the covering count E 'U'x,y ^x,y 

0<x,y<q—l 
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3.8. Implementation and discussion. The MIP approach is easy to im¬ 
plement and also easy to tailor to a search for extreme functions with par¬ 
ticular properties. 

However the approach is limited because floating-point implementations 
are not a good match for finding functions of high arithmetic complexity. If 
q is large, the difference between two slope values Sj often becomes extremely 
small and may completely disappear in floating point fuzz, making the choice 
of the parameter e' difficult. 

Moreover, MIP solvers are generally not the best tool for performing 
an exhaustive search. While listing several solutions should be possible 
by varying the objective function, this is rather difficult to do in practice. 
We have used the solver Gurobi (version 5.6.3) to solve the MIP problem. 
We resorted to setting the Gurobi parameter SolutionLimit=l and calling 
optimize () repeatedly, so that the feasible solutions found by Gurobi before 
reaching the global optimal solution are recorded. 

Despite these limitations, we have obtained new results with the MIP 
approach, which we report in the next two subsections. 


3.9. Result: Optimality of the oversampling factor 3. Using the MIP 
approach described above, we found an example to answer an open question 
in [3j. The open question relates to the following theorem. 


Theorem 3.7 ( [9], Theorem 8.6]). Let m > 3, the oversampling factor, 
be a positive integer. Let tt be a continuous piecewise linear minimal valid 
function for Rf(M./ Z) with breakpoints in -Z and suppose f € ^Z. The 
following are equivalent: 

(1) it is a facet for Rf(R/Z), 

(2) tt is extreme for Rf(R/ Z), 

(3) vr|j_ z is extreme for Rf(^-Z/Z). 

mq ’ 


Dey et al. m gave a function n = drlm_not_extreme_l ( ) that is not 
extreme for i?j(M/Z), but 7r|i z is extreme Rf{- Z/Z). Thus, Theorem 3.7 

does not hold with the oversampling factor m = 1. Basu et al. |7J proved 


Theorem 3.7| for an oversampling factor rn = 4. It was strengthened to 


any m > 3 in [9]. The question whether m > 3 is best possible or can be 
improved to m = 2 was stated in [8, Open Question 8.7]. We resolve this 
question by the following result, which was stated in the introduction as 


Proposition 1.2 


Proposition 3.8. The lower bound m > 3 for the oversampling factor in 


Theorem 3.1 is best possible. Theorem 3.7 does not hold when m = 2. 


Unable to construct an example function by hand, we established this 
result by using the MIP approach. Before we explain the details of the 
search strategy, we describe the structure of the example function that we 
found, kzh_2q_ex ample_l. It is a continuous 4-slope function with q = 37 
and / = see Figure 4. One can verify, simply using the automated 
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extremality test implemented in the software [23], that kzh_2q_example_ 
1 is a non-extreme fnnctioij^j whose restriction to ^Z is extreiiKp^j This 
proves that an oversamplingractor of 3 is optimal. 

In the following, we provide a brief justification of the extremality result 
given by the code. It connects to the theory of equivariant perturbations 
developed in [7], 

Proof of \Proposition 3A . We show that the function n is not extreme for 
Rf(M./7j), by computing the covered intervals. It can be verifiecQ on the 
complex AVi 7 that 


(i) the intervals and [^, indicated by yellow strips in Figure 4 

are uncovered; 

(ii) the others intervals are covered. 

There exists thus a non-zero perturbation function f i, such that ir = \it 1 + 
^7 r 2 where 7T 1 = 7r + ff, n 2 = it — tt are two distinct minimal functions, 
showing the non-extremality of 7r for Rf( M/Z). The function plotted in 
magenta in Figure 4 is such a function 7f. 


By 0 Definition 3.3, Lemma 2.7 and Theorem 4.6], a perturbation func¬ 
tion 7 r is affine linear on the covered intervals, and satisfies 7f(0) = n (/) = 
7r(l) = 0, Ait(x,y) = 0 for any (x,y) such that A-r (x,y) = 0. Consider 
the restriction to ^Z. One can show by linear algebra 24 that the finite di¬ 
mensional linear system has a unique solution 7r| i z = 0. It follows that n 

q 

can only be non-zero on the uncovered intervals, i.e., on [^, and [||, ||]. 
Furthermore, 7r|i z is extreme for R/(- Z/Z). 

Next, we will show that tt \is extreme for Rf{^ Z/Z). Again it can be 
verified that 

(iii) F([^, {^}, [^f, ||]) is an additive edge of AVi z . 

In other words, for [|^, ^], we have A7t(x, t^) = 0, and thus Att(x, = 
0. Following [7], for i e 1, define the translation r t : M —> M, x H > x + t. 
The gold-colored single headed arrow in |Figure 4 indicates the action of 
translation t 4 / 37 , which sends the interval [4y,A] to the interval [ 57 , gf]- 
Therefore, vi(T 4 / 37 (x)) = 7 l(x) + 71 ( 37 ) = n(x) for x € [^ 7 , |=], as 71 ( 37 ) = 0. 
In particular, we have 7 f(||) = 7 f(^|). For r G M, define the reflection 
p r : M — > M, x 1 — > r — x. The gold-colored double-headed arrow in |Figure 4 
indicates the action of reflection pf between the two uncovered intervals, 


21 

22 


As proved by extremality_test (kzh_2q_example_l () ) returning False. 

As proved by simple_finite_dimensional_extremality_test(kzh_2q_example_ 


1(), oversampling=2) returning True. 


23, 


One can type plot_2d_diagram(kzh_2q_example_l () , colorf ul=True) to visualize 
the painting on the complex APi z . 


24 


This could also be verified by simple_finite_dimensional_extremality_test(kzh_ 
2q_example_l() , oversampling=l) which returns True. 
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corresponding to the symmetry condition —x) = 7r(/) = 1 for x E 

[^, We have 7 t(x) + Tt{pf{x)) = 7f(/) = 0 for x G [|y, ^]. In particular, 
7f(|j) + 7r(||) = 0. Therefore, 7f(||) = 7f(||) = 0, the perturbation function 
7 f has values zero at the midpoints of the two uncovered intervals. Since 

•if = 0 on the other intervals which are covered, we have 7f|j_^ = 0- Hence, 

2q 

7r| j_ z is extreme for i?j(^Z/Z). □ 

Our tailored MIP search strategy was to look for functions with properties 
(i-iii) from the above proof. These properties correspond to a prescribed 
partial painting and thus can be expressed by fixing some binary variables. 
We tried out various pairs of ( q , /) for 10 < q < 40. For q = 37 and / = ||, 
we discovered the function kzh_2q_example_lp^| 

3.10. Result: Refutation of the generic 4-slope conjecture. Our 

search also resolves [8] Open Question 2.16] by showing that even for func¬ 
tions whose extremality proof only uses the Interval Lemma, rather than 
the more general techniques from [7] (translations and reflections), many 
slopes are possible. This is in contrast to the first 5-slope functionq 26 | found 
by Hildebrand (2013, unpublished), whose extremality proof requires trans¬ 
lating and reflecting covered intervals. 

Proposition 3.9. There exists a piecewise linear extreme function ir of 
Rj(R/ Z) with more than 4 slopes, such that its additivity domain E(tt) := 
{ ( x,y) : Air(x,y) = 0} is the union of full-dimensional convex sets and the 
lines i 6 Z, i/ £ Z, i+j/G/fZ. 

See the functions kzh_5_slope_fulldim_l etc., which we have made 
available as part of [35] Q They are continuous 5-slope extreme functions 
without any 0-dimensional or 1-dinrensional maximal additive faces except 

^The example can be reproduced using the code in [23] as follows. First, call the 
function write_lpfile_2q(q=37, f=25/37, a=ll/37, kslopes=4, maxstep=2, m=4) to 
generate a MIP problem that maximizes the slope difference S4 — si. The parameter 
a indicates that the uncovered intervals are [a — |,a] and [/ — a, / — a + |], and the 
parameter m decides the heuristic choice of e = t = 1/m. The MIP problem is written 
to the file named mip_q37_f25_all_4slope_2maxstep_m4. Ip. Then, use Gurobi to solve 
the MIP problem and write the solution to the file named solution_2q_example_m4. sol. 
Finally, retrieve the function kzh_2q_example_l from the solution file by calling refind_ 
function_from_lpsolution_2q(’solution_2q_example_m4.sol’, q, f, a). 

2 ®The functions are available in the electronic compendium G5i as hildebrand_5_ 
slope.. . 

^These examples can be reproduced using the code in [23] as follows. First, call write_ 
lpfile(q, f, kslopes, m=12, type_cover=’fulldim’) with appropriate values of q, f 
and kslopes to generate a MIP problem that maximizes the slope difference. For example, 
wesetq=37; f=25/37; kslopes=5. The MIP problem is written to the file named 5slope_ 
q37_f25_fulldim_ml2. Ip. Then, use Gurobi to solve the MIP problem. We set the 
Gurobi parameter SolutionLimit=l and call optimize!) repeatedly, so that the feasible 
solutions found by Gurobi before reaching the global optimal solution are recorded to 
the files named solution_5slope_f ulldim_l. sol, etc. Finally, retrieve the functions 
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Figure 9. The 5-slope extreme function kzh_5_slope_ 
fulldim_l found by our search code (left). Its two- 
dimensional polyhedral complex AV (right), as plotted by 
the command plot_2d_diagram(h,colorful=True) , does 
not have any lower-dimensional maximal additive faces ex¬ 
cept for the symmetry reflection or x = 0 or y = 0 . 


for the symmetry reflections x + y E / + Z and the trivial additivities iGZ, 
y E Z. A graph of the function kzh_5_slope_fulldim_l and a plot of its 
painting on the two-dimensional complex are shown in Figure 9~1 


Remark 3.10. Using our computer-based search we also found extreme 
functions that do have lower-dimensional additive faces, but whose extremal- 
ity proof does not depend on those, 
pies are provided by the functions 
6 _s lop e_fulldim_covers_l etc., which we have made available as part of 

mn 


All intervals are directly covered. Exam- 

kzh_5_slope_fulldim_covers_l, kzh_ 


4. Backtracking search 


4.1. Search via covering paintings. In this section we discuss a new 
search strategy that addresses the limitations of the MIP approach of the 
previous section by using our own implementation of backtracking tailored 
to enumerating covering paintings. 

During our backtracking search, we maintain a prescribed partial paint¬ 
ing as introduced in |Rcmark 3.3 At the root, the minimality conditions 


kzh_5_slope_fulldim_l etc. from the solution files by calling refind_function_from_ 
lpsolution(’solution_5slope_fulldim_l.sol’, q, f). 

no 

" These fulldim_covers examples can be reproduced in the same way as for the 
fulldim examples described in the last footnote, except that type_cover is set to 
’fulldim_covers ’ when generating the MIP problems. For example, the function 
kzh_6_slope_fulldim_covers_l is obtained from the MIP problem 6slope_q25_f8_ 
fulldim_covers_ml2.1p generated by write_lpfile(q=25, f=8/25, kslopes=6, m=12, 
type_cover=’fulldim_covers’). 
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(Theorem 2.2) 

7 To = 0 and 0 < ir x < 1 , 

AlT x ,y = 7T X + Ty — ^{x+y) mod q — 
Tj/ 1 and AiT x ^qf_ x ^ mo( j q 0, 


for x = 1 ,..., q — 1 (Ha) 
Ior0<x<y<q — 1 (lib) 
for x = 0 ,... ,q — 1 (He) 


give the initial prescribed partial painting on AVi 


See 


Figure 10 


for 


an example where q = 4 and / = 4 . To get an extreme function, more 
additivity relations are needed. To achieve this, in our backtracking search 
we successively paint some light grey faces white or dark green, until a 
covering painting is reached. During the painting process, we keep track 
of the consistency of the colors of the faces by using an LP that will be 
explained later in section 4.3 If the current partial painting is infeasible, 


pruning will be performed. 

|Theorem 3.2 and Proposition 2.3 have the following corollary. 


Corollary 4.1. Let it be a continuous piecewise linear function with break¬ 
points in 1Z. If 7 r is extreme for Rf(M./ Z), then there exists a covering 

painting such that 7r|i z is a vertex of the polytope formed by the minimal 
q 

functions for AjQZ/Z) whose additivities correspond to the painting. 

The search for extreme functions ir is thus converted into the search for 
covering paintings. 


4.2. Branching rule. The search tree has a binary structure. Assume that 
we are branching on a triangle F (say the lower triangle whose lower left 
corner is vertex (x,y)) that is colored (light) grey in the partial painting 
of the current node. In terms of|section 31 the value of t xy is undecided. 


Branching on the node creates two children: one where the triangle F is 
painted dark green (£ xy is set to 0) and one where (the interior of) F is 
painted white (£ x ,y is set to 1 ). 

In the child node with dark green F , the following additivity constraints 
hold. 

Att U jV = 0, for every vertex (u, v) of F. (12) 

In the child node with white F, the following strict subadditivity relation 
holds: 

A7 T U v ^ b* 


E 

vertex ( u,v ) of F 


Since the strict inequality constraints are not allowed in linear programming, 
we prefer to replace it by 


^ ' Al T U: y 

vertex ( u,v ) of F 


> £•> 


( 13 ) 


where e is a small positive number. See Remark 3.6 for the value of e. 
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(a) 





implies 

(aa) 


initial conditions: 

0 < 7T* < 1 for x = 1,..., 3, 

7To = 0, 7Ti = 1, 

AtT 0 ,1 — 0, A7T2,3 — 0, 

A7r Xj y > 0 for 0 < x < y < 3. where 

AtTx,!/ ~ "h — 7r ( x+y ) mod 4 




<> 

implies: A7T2,2 = 0, A^i^ = 0 

(bb) 


branch: 1*2,2 
add: A7 t 3 3 


= 0 
<0 



branch: 1*2,2 = 1 
add: A7r 3 3 + 2A7T2 i3 > e 


(<5f 
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4.3. Feasibility and satisfiability checks. To check the feasibility of a 
prescribed partial painting and detect its implied additivity relations, we use 
linear programming. We have investigated two options. As we mentioned 
before, SageMath has a good interface to the Parma Polyhedra Library. 
Its implementation of the double description method provides an attractive 
interface to checking feasibility and for testing implied additivities, all in 
exact arithmetic. This approach, however, appears to be quite slow when the 
dimension of the polytope is large. As a rule of thumb, when the dimension 
exceeds 9, it is better to apply the simplex method instead. Our search code 


uses the GLPK solver, which is integrated well in SageMath (see section A.2 
for details), and allows for warm-starting the simplex method. 

We construct the linear optimization problem as follows. The problem has 
q real variables no ,..., 7 t 9 _i, and some auxiliary variables A 7 x x ^ y (0 < x < 
y < q — 1) that represent the subadditivity slacks. The initial constraints on 
the variables are given in & The subadditivity and additivity constraints 
are reflected by the bounds of their slack variables^] These bounds will vary 
along the backtracking process. If we walk downwards in the tree, then we 
append either (12) or (13) to the constraint system; conversely, if we walk 
upwards in the tree, then we remove the constraint. 

Such changes in the constraint system could affect the feasibility of the 
problem. Due to the update of the variable bounds, the dual simplex method 
starting off the last basis is called to check whether the partial painting 
remains feasible at the current node. If the node is infeasible, the whole 
sub-tree will be pruned. For example, the node (f) in Figure 10 is infeasible, 
as < 0 and A 773 33 > e > 0. 

The current constraint system may imply some new additive faces in the 
prescribed partial painting. To check if a vertex ( u , v) is implied additive, 
we call the primal simplex method starting off the last basis to maximize 
the objective function Air u , v . If the optimal value is 0, then A 7 r U) „ = 0 and 
thus the vertex ( u , v ) is implied additive (dark green). In terms of section 3 
Pu,v = 0. Subject to the inclusion constraints ([ 6 ]), the colors of edges and 


triangles are updated accordingly. See the nodes (bb) and (dd) in Figure 10 
for examples. 


4.4. Heuristic choice of the branching triangle. By the invariance of 

the subadditivity condition under exchanging x and y , only the upper left 

triangular part of the complex AThg needs to be considered for painting. 

q 

In our experiments we found that an exhaustive search of all paintings 
of the two-dimensional complex is too expensive even for a moderate size 
of q. Thus, to reach a covering painting quickly, we branch on the colors of 
the triangles (£ x>y and u X)V ) rather than on the colors of the vertices (p x , y )- 
Furthermore, a heuristic painting strategy is applied: While picking a light 
grey triangle to paint dark green or white in branching, we consider those 

nq 

We need to introduce these slack variables explicitly due to limitations of warm¬ 
starting in the SageMath interface. 
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triangles F whose projections pi(F) and p 2 (F ) are currently uncovered. 
This is of course restrictive, and so we cannot guarantee that our search 
code will find all covering paintings and all extreme functions. However, it 
has proved to be a successful heuristic strategy. 

At each node, we choose one candidate triangle F (i.e., one i x ^ y or u x>y 
variable) among all these considered triangles. It is defined as the smallest 
triangle in lexicographical order, whose color has not been branched on yet 
in the search tree. 


4.5. Incremental computation. For the purpose of improving the run¬ 
ning time efficiency, all computations in the backtracking search, such as 
updating covered intervals, are done in an incremental manner. 

More precisely, we maintain a list of connected components that we have 
mentioned in section 3.1 When a new triangle F is painted as additive in 


the prescribed partial painting, the components that contain the projection 
pi(F) or p 2 (F ) or ps(F) are merged into one big component, and all intervals 
in this new component become covered. When a new edge (one-dimensional 
face) F is painted as additive (dark green), the components that contain 
its projection intervals are merged into one big component. If the new 
component contains an interval that was covered, then all intervals in this 
new component are covered. In such a way, the new covered intervals after 
adding an additive face can be computed incrementally from the covered 


intervals in the previous step. For example, the node (a) in Figure 10 has 
three connected components C\ = {[0, |]}, C 2 = {[ 5 , 5 ], [|, 1]} and U = 
{[2 ’ 3 ]}’ w h ere the first two are covered and the last one is not. When 
£ 2,2 is set to 0 in its child (b), U is merged into C 1 , and becomes covered. 
Thus, the node (b) has two connected components C[ = {[0,|],[|,|]} and 
C2 = {[4, 5], [f , 1 ]} that are both covered. 

We mentioned above that a function ir whose additivities satisfy the pre¬ 
scribed partial painting has the same slopes on the intervals in each compo¬ 
nent, see fTj Remark 3.6]. By counting the connected components, we get an 
upper bound on the number of slopes that the function n could have. This 
allows us to prune subtrees that cannot contain functions with the desired 
number of slopes. 

The knowledge of connected components is also used at the end of 


ver¬ 


tex filtering mode” (section 2), to check efficiently whether all intervals are 
covered, as follows. If an interval [|, does not belong to any covered 
connected component from the partial painting and if none of the values Q 
is 0, then 7 r has uncovered interval and hence is not extreme. This saves 
us from running the full extremality test on n, which would consume more 
time according to Figure 5]-middle. 


4.6. Heuristic search algorithm. We summarize the backtracking search 
via covering paintings in Algorithm 2 
mode” search in our code. 


It is referred to as the “heuristic 
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Algorithm 2: heuristic mode 


( 1 ) 

(2) 

( 3 ) 

( 4 ) 


( 5 ) 

( 6 ) 


The root node is the initial painting given by the minimality 
conditions ©; 

Decide for a candidate triangle F of the painting using covered 
intervals; 

A node is branched into two child nodes: 

For the child node in which F is additive (dark green), 

• add the new additivity relations © to constraints; 

• look for implied additive vertices and triangles; 

• update covered intervals; 

• if the node is infeasible, backtrack; 

• if the number of the connected components is less than the 
desired number of slopes, backtrack; 

• if a covering painting is found, output it and backtrack; 
For the sub-node in which F is strictly subadditive (white), 

• add the strict subadditivity relation (13) to constraints; 
Traverse the search tree in depth-first order. 


For each covering painting returned by the algorithm, we use the mini¬ 
mality conditions © and the additivity relations Att^ = 0 specified by 
the green vertices (x, y) in the painting, to construct a polytope. (The 
strict subadditivity constraints Air x<y > 0 corresponding to the green ver¬ 
tices (x,y) in the painting are neglected.) We enumerate the vertices of the 
polytope. By Remark 3.5| and |Theorem 3.2[ the interpolation 7r of a vertex 
is extreme for Rf(M./ Z). 


7T| 1 

q. 


4.7. Combined search algorithm. The above search algorithms work 
well for relatively small </, but become inefficient when q is large: The vertex 


filtering search (Algorithm 1) wastes time on enumerating numerous vertex- 


functions in high dimension, most of which are non-extreme for the infinite 


group problem; the heuristic backtracking search via covering paintings (Al¬ 


gorithm 2) suffers from the combinatorial explosion in branching and the 


general performance penalty from using Python. 

We propose to combine the vertex filtering search and the heuristic back¬ 
tracking search to obtain a better performance. The combined algorithm 
starts with branching, but outputs the partial painting and backtracks at a 
certain depth before reaching a covering painting. For each generated partial 


painting, the algorithm constructs the polytope as described in section 4.6 


and then performs a vertex enumeration as described in the vertex filtering 
search. It remains to determine a good stopping criterion for branching. 

Since the vertex enumeration algorithm has a good performance for low¬ 
dimensional polytopes, we wish to use the dimension as the stopping cri¬ 
terion. However the actual dimension of the polytope given by a painting 
is unknown, unless it has been constructed, when it is too late. Therefore, 
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instead of the actual dimension, we use expected dimension as the stopping 
criterion in our code. This expected dimension can be computed efficiently 
without calling PPL to construct the polytope. We set up a g-column ma¬ 
trix (cs_matrix) to record the equality constraints on (wo, ■.., n q ), which 
are ttq = 0, the symmetry constraints and the additivity constraints spec¬ 
ified by the painting. The matrix is maintained dynamically during the 
backtracking process. Define the expected dimension to be the co-rank of 
the equation system: exp_dim := q — rk(cs_matrix). 

The algorithm switches from backtracking to vertex enumeration once 
the expected dimension becomes smaller than a certain threshold. Table 4 


shows that a value around 11 is the best empirical threshold for finding an 
extreme function with many slopes quickly. 

The combination of vertex filtering search and heuristic backtracking 
search produces a more powerful search algorithm, which is called “com¬ 
bined mode” in our code. We summarize it as Algorithm 3 


Algorithm 3: combined mode 


(1) Run the heuristic search (Algorithm 2), with one more stopping 
criterion added to its step [4j 

• append the new equations to cs_matrix; 

• compute exp_dim := q — rk(cs_matrix); 

• if exp_dim < threshold (empirically, threshold = 11), 
output the partial painting and backtrack; 

(2) For each partial painting returned by phase [T] construct the 
corresponding polytope and run the vertex filtering search 


(Algorithm 1, steps 3-4). 


4.8. Results. Using the combined 
able to find up to 7-slope extreme 

slope_l to kzh_7_slope_4. 


search (Algorithm 3), 
functions for q < 34, 


our cod^] was 
namely kzh_7_ 


Of) 

ou By running search_kslope_example (k_slopes, q, f, mode=’combined’) with 
various values of k_slopes, q, f. For example, kzh_7_slope_l can be obtained by set¬ 
ting k_slopes=7; q=33; f=ll. 












Table 4. Vertex enumeration in high dimension vs. Combinatorial explosion in branching 
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5. Targeted search for extreme functions with many slopes 


We observed that many of these newly discovered extreme functions with 
many slopes possess the invariance: / = 1/2 and 7 Tj = n q -i (0 < i < q/ 2 ). 
Targeting the search to functions with the invariance property allowed us to 
find more new extreme functions with many slopes whose values of q were 
twice as large as before. In addition, their painting on the complex AVij 


often includes special patterns, as shown by Figure 11 for example. 


5.1. Construction of a family of prescribed partial paintings. We 

then targeted the search to functions for larger values of q with prescribed 

partial paintings that mimic these patterns. Let q = 36r + 22, where r E 

Z, r > 1. We construct the prescribed partial painting on the complex APi z 

q 

in 2 (r + 2 ) steps as follows. 

Step 0: Paint the lower triangles whose lower left corners are the ver¬ 
tices i 0 ! ( ° 1 ( ° 1 ('(<J-2)/4g' ) ((q-2)/4qs , ((q-2)/2q\_ , . . 

tlCeS l 0 ) ’ l (q - 2 )/ 4<7 / ; l (q -2)/2q ) ’ l (q-2)/4q ) ' V 0 Jand^ q ) j tills IS 

illustrated in Figure V2\ -l. 

Step 1: Paint the upper triangles whose upper right corners are the ver¬ 
tices (*y q q ) for (®.) = ( 9r ^_ 5 ), (^+ 5 ) and ( 9r ' 2 h5 ). Paint the lower triangles 

whose lower left corners are the vertices ( l / q q ) for (*) = ( 9 /^ 4 ), ( 9r + 4 ), 

/ 4 \ / 6 \ / 9r+2 \ / 9r+4 \ / 9r+2 \ / 9r+4 \ / 9r+4 \ / 9r+2 \ , / 9r+2 \ _ 

l 9r +2 ) ’ V 9r+2 / 1 V 9r+4 ) ’ V 9r+4 ) ’ V 9r+2 ) ’ V 9r+2 ) ’ V 4 J’V 4 ) anQ l 6 / ’ 

see |Figurc 12| -2 . 

Steps t = 2, 3,..., r: Paint the parallelogram whose vertices are () for 


• 6t—9 

. j / v 9r —3t+10 
shown in 


;) = 


Figure 12 

9r—3i+5 


) , ( 9 r—3t+10 ) > ( 9r—3t+5 ) ’ ( 9r-3f+5 ) With the 0ran § e P attem 
3(a)). Paint the square whose vertices are (®.^) for 


/z\ / \)r—6t-\-b \ / yr—cSt+lU \ (vr—6t+b\ (\)r-dt-\- 1U \ -,1 ,i , 

(j) = ( 9i—3t +10 ) ’ ( 9r—3 f+10 ) ’ ( 9r—3t+5 ) ’ ( 9i-3t+5 ) Wlth the ° ran S e P at ~ 
tern shown in|Figure 12-3(b). Paint the parallelogram whose vertices are 

($) fa (’) =Tli l? 5 ). (V -T). (^I 10 ), CuT) with the or- 
ange pattern shown in Figure 12]-3(c). 


Step (r +1): Paint the triangle whose vertices are( ) for (\) = ( 
( 6 r+ 7 ) ’ (6r+3 ) w ^h the red pattern shown i 


Figure 12-4 


6r—3 \ 
6r+7 ) ’ 


Step 0 to Step (r + 1) construct the painting on the lower left triangular 

part of the complex AVi%. The painting on the upper right triangular 

q 

part of the complex APi z will then be determined through a mapping. 
Specifically, in Step t for t = r + 2, r + 3,..., 2r + 3, we paint the triangles 
whose images under the mapping (^) i —> (j_^) are colored in Step (2 r + 

3 — t). We also paint the diagonal lines {( *): x + y = \,0 < x < \} and 
{(y) ■ x + y = ^,^ < x < 1}, which correspond to the symmetry condition 
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Figure 11. Special patterns on the two-dimensional poly¬ 
hedral complex AVi z . Left, the AVi z of the 6 -slope ex- 
treme function kzh_6_slope_l with q = 58. We observe 
that the additive triangles are located in the lower left and 
upper right corners. The function has the same slopes on 
the intervals that are projections of the same color additive 
triangles. The 6 -pointed star patterns appear several times. 
Right, the lower-left corner of AVi% of the 10-slope extreme 


function kzh_10_slope_l with q = 166, where we see that 
the 6 -pointed stars are actually the result of additivity pat¬ 
terns within certain intersecting quadrilaterals {black), which 
connect like links of three chains. The detailed structure is 


described in Figure 12 


of minimal valid functions. In the following, we shall refer to the painting 
constructed as above on AVi z as the prescribed partial painting. 

5.2. Properties of functions satisfying the prescribed partial paint¬ 
ings. Recall the notion of connected component discussed in |section 4.5 
Connected components of a painting on AVi z are disjoint subsets of {[^, ^-] 

i = 0,1,... ,q — 1}. They satisfy the following properties. If F is an ele¬ 
mentary upper or lower triangle whose vertices are colored on the paint¬ 
ing, then pi(F),p 2 (F) and ps(F) are in the same connected component. 
If F is an elementary horizontal, vertical or diagonal edge whose vertices 
are colored on the painting, then the two intervals among its projections 
pi{F),p 2 {F),ps(F) are in the same connected component. In particular, the 
colored diagonal lines corresponding to the symmetry condition yield that 
l\' an< ^ [</’ are * n the same connected component, where i + j = 
§ - 1 = 18r + 10. 
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Pi 


P 2 

18r+10 1 2 

5 "25“ 


r 


1/2 


9r+5 _g-2 b... 


-u— 




-//- 




(9-2)/4? (q~2)/2q 

P3 


AT 70 of Step 0 


Pi 



(' St-9)/q 


(9r-34+10)/g 


P2 

(9r—34 +9)/g 
(9i—34 +7)/q 
(9r—34 +5)/g 



AT 7 * of Step t 

Pi 


P3 

(18r—6i +19)/g 

(18r—6f +17)/g 

(18r-64+15)/g 

(18r-64+14)/g 

(18r-64+12)/g 

(18r—64 +10)/g 

(9r+34+5)/g 
(9r+34 +3)/g 
(9r + 34+l)/g 


(Gr—S)/q (q—4)/6q (6r + 7 )/g 



P3 

(12r+13)/g 

(12r+ll)/g 

(12r+9)/g 

(12r+8)/g 

(12r+7)/g 

(12r+6)/g 

(12r+4)/g 


AP r+1 of Step (r + 1) 


Figure 12. Prescribed partial paintings. 
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Lemma 5.1. Let q = 36r + 22, where r E Z, r > 1. Let the prescribed partial 

painting on AVi z be constructed as above. Then all intervals are directly 
<? 

covered. The prescribed partial painting induces exactly 2(r + 2) connected 
components. 


The proof of this and the following results appears in Appendix C 
1, q = 36r + 22 and / = 1/2 


Suppose r E Z, r > 1, q = 36r + 22 and / = 1/2. Let II r be the set 
of continuous piecewise linear minimal valid functions n with breakpoints 
in l q Z, satisfying in addition the invariance condition tt(x) = 7 r(l — x) for 

0 < x < | and A 7 t(x, y) = 0 for any 0 < x, y < 1 such that the point (x , y ) 
is colored in the prescribed partial painting on AVi z . 


Clearly we can describe the functions 7r E II r in the space of the vari¬ 
ables (7To,7Ti,... ,TT q )] let us denote the corresponding polytope by V r . Due 
to |Lemma 5.1 [ we can describe them also in the space of the slope values 
(so, si,..., s r _|_i) on the connected components; let us denote the correspond¬ 
ing polytope by S r . (Note that st = —S 2 r+ 3 -t by the invariance condition, 
for t = r + 2, r + 3,..., 2r + 3.) 


Lemma 5.2. There is a linear isomorphism between the polytope V r in the 
(-To, 7Ti,..., 7 T q ) variables and the polytope S r in the (so, si,..., s r +i) vari¬ 
ables. 


The proof (in Appendix C) gives an explicit mapping between the values 


7 Tj and the slopes s*. With the help of these formulas, we can then show 
that the slope values are non-increasing. 


Lemma 5.3. Let so, si,..., s r+ i be as above. Then so > s\ > ■ ■ ■ > s r+ 1 . 

The computer-based search can now be run in the space of the slope 
variables (so, si,..., s r +i), which has the benefit of having a much lower 
dimension. However, if r is large, the search is still nontrivial. In order 
to speed up the search we prescribed extra additivity constraints. This 
approach was successful in finding extreme functions with up to 28 slopesp*] 

Our search also revealed that in general, we cannot expect the existence of 
an extreme point for which the sequence of slope values is strictly decreasing 
(so > si > ■ ■ ■ > s r -|_i). However, we conclude this section with a weaker 
conjecture. 

Conjecture 5.4. There exists an extreme point of the polytope S r with H(r) 
different slope values s*. 


■^By running pattern_extreme(r, k_slopes) with various values of r and k_slopes. 
For example, kzh_28_slope_l can be obtained by setting r=21; k_slopes=28. In order 
to reduce the running time of the targeted search when r is large, the code pattern, 
extreme () imposes some extra colored vertices on the prescribed partial painting. They 
correspond to extra additivity constraints on the function n, which are often satisfied by 
the previously discovered many-slope extreme functions. Concretely, when r > 16, we 
assume that An x , y = 0 for (x, y) = (6 r + 5, 36r + 18), (6r + 7, 36r + 10), (6r + 7, 36r + 
12), (6r + 10, 36r + 3), (6r + 11, 36r), (9r - 18, 9r - 18), (9r - 12,9r - 12), (9r - 9,9r - 
9), (9r - 3, 9r - 3), (9r + 3,9r + 3). 
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We abandoned work on this conjecture in late July 2015 when Basu et al. 
[6] announced a different construction that gives extreme functions with an 
arbitrary prescribed number of slopes. 

5.3. Result: Extreme functions with many slopes. The targeted search 
was very successful in finding functions with large numbers of slope^} We 
thus obtained the following result, which we have stated already in the in¬ 
troduction. 

Theorem 5.5. There exist continuous piecewise linear extreme functions 
with 2, 3, 4, 5, 6, 7, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, and 28 slopes. 


Appendix A. Implementation details 


In this appendix, we describe some aspects of our implementation in Sage- 
Math [34], an open-source mathematics software system that uses Python 
and Cython as its primary programming languages and interfaces with var¬ 
ious existing packages. 


A.l. SageMath interface for vertex enumeration. PPL is a standard 
package in SageMath and comes with an efficient Cython interface. Our code 
uses in particular the class C_Polyhedron for computations with closed con¬ 
vex polyhedra. A polytope can be built starting from a system of constraints 
cs of class Constraint_System via p = C_Polyhedron(cs) , where the con¬ 
straint system cs is a finite set of linear equality or inequality constraints 
(class Constraint). One calls p.minimized_generators () to enumerate 
the vertices of p. PPL also allows for feasibility checks and satisfiability 
checks (see section 4). The feasibility check can be realized by calling p. is_ 
empty (). The satisfiablity check efficiently tests whether a given inequality 
or equation c is satisfied by all points in a polytope p. It is accessed by calling 
p.relation_with(c).implies(Poly_Con_Relation.is_included!)). 

Once lrslib EU2] has been installed as an optional package in SageMath, 
it is possible to call the programs lrs and redund from SageMath. Our code 
includes a SageMath interface that reads or writes polytopes in the lrslib 
format. The lrslib command redund can thus be used in conjunction with 
PPL as a preprocessor for vertex enumeration. 


A.2. SageMath interface for linear programming. We pointed out in 
the necessity of using an LP solver with warm-starting capabil- 
feasibility and satisfiability checks described above in the high¬ 
dimensional case. We use the SageMath class MixedlntegerLinearProgram 
as an LP modeling system. Within this framework, a new LP problem m 
can be created by m = MixedlntegerLinearProgram(maximization=True, 
solver = "GLPK"), requesting the GLPK solver as its numerical backend. 
In contrast to other backend implementations, including PPL’s rational LP 

oo , , 

We have made the functions available as part of the Electronic Compendium [35] as 
kzh_7_slope. . . and kzh_28_slope. . ., etc. 


section 4.3 
ity for the 
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solver, the GLPK backend has the crucial warm starting capability. We call 
v = m.new_variable(real=True, nonnegative=True) to define a Python 
dictionary v of non-negative continuous variables for the problem m. The 
upper and lower bound of a variable, say v[ 0 ], can be changed via m. set_ 
max(v[ 0 ] , max) and m. set_min(v [ 0 ] , min) respectively. If the variable 
is unbounded above or below, then one sets max=None or min=None respec¬ 
tively. The method m.add_constraint(linear_function, max, min) 
sets up a new constraint min < linear_function < max for the problem m. 
The objective function of m is defined by m. set_obj ective (obj). For a fea¬ 
sibility check, we can use obj=None. We request that GLPK use the simplex 
method to solve the LP via m. solver_parameter (backend.glp_simplex_ 
or_intopt, backend.glp_simplex_only). According to the setting 
m.solver_parameter("primal_v_dual", "GLP_PRIMAL") or m.solver_ 
parameter("primal_v_dual", "GLP_DUAL"), the primal or dual simplex 
method is applied respectively. We call m. solve (objective_only=True) 
to solve for the optimal value. If it signals a MIPSolverException, then the 
problem is infeasible. 

References: 


Appendix B. Limitations of search based on q x v grid 

DISCRETIZATION 


In this section, we discuss limitations of the search based on q x v grid 
discretization, an alternative search strategy that was used by Chen (121 and 
Hildebrand (2013, unpublished). 

Consider continuous piecewise linear functions n: M/Z — > [0,1], with 
breakpoints in ^Z for some q G Z + and 7 r( 0 ) = 0. Suppose without loss 
of generality that / G ^Z. 

As mentioned in |section 1 .4] there are two natural ways to discretize the 
space of functions 7 r: discretizing function values 7 Tj = 7 r(|) for i G {0,..., g} 


Figure 2 


G {1, ..., q}. See again 
The following lemma shows that they are equivalent. 


and discretizing slope values qsi on for i 


Lemma B.l. Let v be a positive integer. The following are equivalent: 

(1) 7 Tj G \7L for each i G {0, ..., q}. 

(2) Si G \7L for each i G {1,... ,q}. 

Proof. Since ttq = n q = 0 and Si = 7r* — 7Tj_i for i = 1 ,... ,q, the lemma 
follows. □ 


B.l. A lower bound on (a proxy for) arithmetic complexity. In the 

following, we investigate the worst-case complexity of the search based on 
qxv grid discretization, by estimating the arithmetic complexity of extreme 
functions, i.e., the largest value v needed for any extreme function ir with 
breakpoints in ^Z. 
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We use the notations d ex t, d ve r , and dbas, satisfying d ex t < d ve r < dbas 5 
and the constraint system Ax. = b of the polytope IIj(^Z/Z), which were 
introduced in section 2.6 It is hard to determine the precise value of d ex t as 
a function of q, and we are not able to show an exponential lower bound for 
it. For estimating the growth rate of d ex t, we are satisfied with a simplified 
study, using db as as a proxy. We show the following exponential lower bound. 

Lemma B.2. Let q > 3 be an odd positive integer. Let f G 0 < / < 
1, such that qf and q are coprime integers. Let Ax = b, x > 0 be the 
constraint system of Theorem 2.2 written in the standard form. Then d^ as , 
the maximum absolute value of the determinants of simplex basis matrices 

q — 1 

of A, is at least 2 2 . 


Proof. It suffices to show the existence of a basis matrix B of A with 

q — 1 

|det-B| > 2^~. To find such a B, we first prove the following claim. Be¬ 
cause q is odd, the operation of multiplying by 2 (mod q) is invertible. For 
x G {0,1,..., q — 1}, denote the unique y G {0,1,..., q — 1} satisfying 2 y = x 
(mod q) by x/2. 

Claim B.3.Q Let q and f be as above. There exists a sequence (ao, ai,..., 
a q - 1 ) of integers with ao = 0 , a± = qf and 02 = qf/2 (mod q) such that the 
following conditions hold: 

(1) for odd i > l we have a* = ajf 2 (mod q ) for some j < i; 

(2) for even i > 2 we have ai = qf — a*_ 1 (mod q). 

(3) {a 0 ,ai,... ,a ? _ 1 } = {0,1,..., q - 1}. 

Proof. We construct the sequence as follows. Suppose that ao, ai,..., ak 
are determined for some even k > 2 , such that conditions ( 1 ) and ( 2 ) are 
both satisfied for i < k, and that ao, at, ■ • ■, a*, are all distinct. Let S = 
{ao, oq,..., a*;}. We choose ak +1 and ak +2 by selecting an element s G 
S such that s /2 (mod q) ^ S, and then taking a^+i = s /2 (mod q) and 
ak +2 = Qf — s/2 (mod q). It suffices to show the existence of such an s G S 
whenever S 7 ^ Z/qZ. 

Suppose that s/2 G S for every s G S. Since qf G S and qf and q 
are coprime, S must contain the coset qfH = {qfh : h G H}, where H 
is the multiplicative subgroup of (Z/gZ)* generated by 2. In particular, 
qf, 2qf G S. 

By conditions ( 2 ) and then ( 1 ), we deduce that S also contains qf — qfH 
and (qf — qfH)H = qfH — qfH. Applying this argument repeatedly, we see 
that S contains qfH — qfH + qfH — ■ ■ ■ ± qfH for any number of iterations. 
Since 1,2 G H, qfH contains qf and 2 qf, and thus any multiple of qf can 
be written in the form qfH — qfH + qfH — • • ■ ± qfH. Since qf and q are 
coprime, we conclude that S contains all of Z/qZ. □ 


Define the row vectors Ro, R\,..., R q -\ G Z 9 using the sequence a 0 , 01 ,..., 
a q - 1 constructed above, as follows. Let Ro be the row vector with the only 

■^Thanks go to Xuancheng Shao for the help in proving this claim. 
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nonzero entry 1 appearing in the column indexed by ao = 0 , corresponding 
to the constraint 7 To = 0 . Let Ro be the row vector with the only nonzero 
entry 2 appearing in the column indexed by 02 , corresponding to the sym¬ 
metric constraint 7 r a2 + 7 r a2 = 1. For i = 1 and i = 4, 6 ,..., q — 1, let Ri be 
the row vector with two nonzero entries 1 appearing in the columns indexed 
by a,; and a,_ 1 , corresponding to the symmetric constraint n ai + 7 T ai _ 1 = 1 . 
Finally, for i = 3, 5, ..., q — 2, let Ri be the row vector with nonzero en¬ 
try —2 at index a* and entry 1 at index 2 a* (mod q), corresponding to the 
subadditive constraint 7 r ai 4 - 7 r ai > 7 T 2ai mo d q- 

The basis matrix B is obtained by taking the slack variables for the subad¬ 
ditivity constraints R 3 , R §,..., R q -2 in A as non-basic variables and others 
as basic variables. To compute det B , first expand out the columns corre¬ 
sponding to slack variables. We are left with a qx q matrix B' consisting of 
the rows Rq, R\, ..., R q ~i, and |det B\ = |det B'\. See 
case of q = 11, / = 3/11. 

To compute detiT, start by expanding along the row Rq containing a 
unique nonzero entry 1 and end up with a new matrix with this row and 
column a 0 removed. In the second step, expand along R\. noting that the 
only nonzero entry remaining in this row is 1 at column a\. We arrive at a 
new matrix with this row and column a\ removed. In the third step, expand 
along i? 2 , noting that the only nonzero entry remaining in this row is 2 at 
column a 2 . We then arrive at a new matrix with this row and column a 2 
removed. In general, during the (k + l)-st step, we expand along the row Rk 
which contains a unique nonzero entry at column Ofc, whose value is either 
—2 or 1 depending on whether k > 3 is even or odd. 

The computation terminates in q steps. It follows that det B' is equal to 
the product of all these unique nonzero entries, (q — l )/2 of which are ±2 

q — 1 

and the remaining are 1. Thus |deti?| = 2 ~. □ 


Example B.4 for the 


Example B.4. Consider the case q = 11,/ 
the sequence 


11 ■ 


By 


Claim B.3 


we have 


(ao, ai,..., aio) 


(0,3,7,9,5,10,4,8,6,2,1). 


The following matrix B' is a q x q submatrix of A, where the row Ri corre¬ 
sponds to the: 

• constraint tto = 0, for i = 0; 

• symmetric constraint 7 r a2 + ir a2 = 1 , for i = 2 ; 

• symmetric constraint 7 r ai + 7 r ai _j = 1, for i = 1,4, 6 ,... , q — 1; 

• subadditive constraint — 27r Qi + 7 r 2 ai mod g < 0, for * = 3, 5, — 2 

°f n /(JZ/Z). 
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B' = 


01 234567 8 

/ 1 0 000000 0 

1 0 0 1 0 0 0 0 0 

00 000002 0 

00 000001 0 

00 000100 0 

00 000000 0 

00 001000 0 

00 000100 -2 

0 0 0 0 0 0 1 0 1 

00 -2 01000 0 

V 0 1 100000 0 


9 

0 
0 
0 
-2 
1 

1 -2 


10 

°\ 

0 

0 

0 

0 


0 

0 

0 

0 

0 


1 
0 
0 
0 
0 ) 


Ro 
Ri 
Ri 
R3 
Ra 
Re 
Re 
R7 
Rs 
Rg 
Rio 


B.2. An upper bound. We now prove the exponential upper bound stated 
in Isection 2.61 


Proof of Lemma 2.6. To compute detB, we first expand out the columns 
corresponding to slack variables as in the proof of Lemma B.2 We are 
left with an n x n matrix B', where n < q. Denote the rows of B' by 
Ro, R \,..., R n —i G IP. We distinguish the types of constraints that the 
rows Ro, ili,, Rn-i £ correspond to in Theorem 2.2 
constraint 7 To = 0 , for which 
Tr T ^(qf—x) mod q 

of n and qf, for which ||i?j ||2 < V 2 if x / (qf — x) mod q and ||i?i ||2 < 2 


There is one 

| Ri 11 2 = 1; there are m symmetric constraints 
= 1 , where m = |, or depending on the parities 
12 < ^ if x / (qf — x) mod q and 
if x = (qf — x ) mod q ; and there are n — 1 — m subadditive constraints, for 
which 11 Ri 11 2 < y/5. Therefore 


n—1 


|det B\ = |det B'\ < < (Vb) 71 -^ 1 (V2) m+2 < 10 9/4 . 


□ 


i=0 


B.3. Conclusion. Although the question is not conclusively settled, Lem¬ 


mas 


B.2| and |2.6| indicate that the value v needed in the q X v grid dis¬ 
cretization grows exponentially with q. The empirical results of <i e xt and 
r obtained by the vertex filtering search (see section 2 ) confirm this ex¬ 


ponential growth, as shown in |Table~3| and Figure 6 


We conclude that the search based on the q x v grid discretization for 
breakpoints and function values (or, for breakpoints and slope values, by 


Lemma B.l) is not suitable for an exhaustive search if q is large, due to its 


high worst-case complexity. 


Appendix C. Proofs of the theorems in Isection ~5l 


Proof of Lemma 5.1 
step, for t = 0,1,..., 2r + 3. 


Let AV t denote the set of colored triangles in the f-th 
Consider their projections pi(AV t ), p 2 (AV t ) 
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and p 3 (A'P t ). Let p(AV 1 ) := U£=i PkiAV*). By Figure 12 
Pi(AV°) = p 2 (AV°) = p 3 (AV°) = {[|, : * = 0, 9r + 5,18r + 10}; 

Pi(AV 1 ) = p 2 {AV l ) = {[i,if]:i = 1 , 2,4, 6 , 9r + 2, 9r + 4}, 
p 3 (AV l ) = {[J, : i = 9 r + 6 , 9r + 8,18r + 4,18r + 6 ,18r + 8,18r + 9}; 

PiiAV*) = p 2 (AV t ) = {[i,i±^] : i = 6 t — 9, 6 t — 7, 6 f — 5, 6 t — 4, 6 t — 2 , 6 f, 
9r — 3t + 5, 9r — 3t + 7, 9r — 3f + 9}, 

p 3 (A'P t ) = {[|, i±i] : i = 9r + 3t + 1 ,9r + St + 3,9r + St + 5,18r -6 1 + 10 , 
18r -6 1 + 12,18r -6 1 + 14,18r -6 1 + 15,18r -6 1 + 17,18r -6 1 + 19}, 
for t = 2 ,3,..., r; 


pi(AV r+1 ) = p 2 {AV r+1 ) = {[^, : i = 6 r - 3,6r - 1, 6 r + 1, 

6 r + 2, 6 r + 3, 6 r + 4, 6 r + 6 }, 
p 3 (AV r+1 ) = {[J, *±i] : i = 12r + 4,12r + 6 ,12r + 7, 

12r + 8,12r + 9,12r + 11,12r + 13}. 


The sets p(A'P t ) for t = r + 2 , r + 3,..., 2r + 3 can be obtained through the 
mapping i o 1 - i. 

Note that for each t = 0,1,..., r + 1, [|, G pi(A'P t ) = p 2 (A'P t ) if 
and only if [^ - ^, i - J] = [ 18r + 10 ^ , 18r + 11 ~ i ] G p 3 (A7 ?t ). Therefore, 
the set p(AV t ) is stable under the reflection corresponding to the symmetry 
condition. 

We now show that, for each t = 0,1,..., r + 1, the intervals in p{A'P t ) 
are from the same connected component. 

In Step 0, consider the 3 green triangles on the yellow diagonal stripe with 
p 3 = Since they have the same p 3 projection, their p 2 projections 

which form the set p(AV°) are from the same connected component. 

In Step 1, consider the 6 green triangles on the yellow diagonal stripe 
with p 3 = [^± 6 , ^±1], Their p 2 projections V := {[^, ^±1] : i = 1,2,4,9r + 
2,9r + 4} are from the same connected component, say C. Let F denote 
the green lower triangle whose pi(F) = [ 9r + 2 , and p 2 (F) = [^, 

Then p 2 (F) G C since pi(F) ePC C. p±( AV 1 ) = p 2 (A'P 1 ) C C. Using the 
reflection x ha — x) mod 1 corresponding to the symmetry condition, we 
have p( AV 1 ) C C. 

In Step t for t = 2,3,..., r, consider the 8 orange triangles on the yellow 
diagonal stripe with p 3 = [ 9r + 8t + 3 , 9r+3t+4 ]- Since they have the same p 3 
projection, their p 2 projections V := {[^, : i = 6 t — 7, 6 t — 5, 6 f — 4, 

6 1 — 2, 9r — 3t + 5, 9r — 3t + 7, 9r — 3t + 9} are from the same connected 
component, say C. Let F\ denote the orange upper triangle whose p\ (F \) = 
[ 9 r ~f +9 , 9r ~ 3 g t+1 ° ] and p 2 (Fi) = [%=*, ^]. Let F 2 denote the orange lower 
triangle whose Pl {F 2 ) = [ 9r ~ 3t+5 , 9r ~ 3f+6 ] and p 2 (F 2 ) = Since 
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pi(F\) e V C C, p 2 {F\) € C. Similarly, since p\(F 2 ) ePC C, p 2 (F 2 ) £ C. 
Thus pi (AT 3 *) = p 2 (AV t ) C C. Using the reflection x i-A (| — x) mod 1 
corresponding to the symmetry condition, we have p(A'P t ) C C. 

In Step r + 1, consider the 4 red triangles on the yellow diagonal stripe 
with P 3 = [ 12 ^ +4 , 12 g +5 ]- Their P 2 projections R := {[g> : i = 6 r — 

3, 6 r — l, 6 r + 4, 6 r + 6 } are from the same connected component, say C. 
Let F k denote the red upper triangle whose pi(F k ) = [ br + k , 6r +^'+ 1 ] and 
P2(F k ) = [^,^±5], for k = 1,2,3. Since p 2 (F k ) e V C C, pi(F fc ) € C 
for each k = 1,2,3. Thus pi(AV r+1 ) = p 2 (AV r+1 ) C C. Finally, by the 
reflection x 1 —> — x) mod 1 corresponding to the symmetry condition, 
we conclude that all elements of p(AV r+1 ) are from the same connected 
component C. 

By construction, p(AV°), p(AV 1 ),... ,p(AV r+1 ) form a partition of the 

set {[|, Al] : i = 0,1_ 18r + 10}. Recall that q = 36r + 22. Then by the 

invariance under the mapping x i-a 1 — x, we have that p(AV r+2 ), p(AV r+5 ), 
... ,p( AV 2r+3 ) form a partition of the set { [|, A] : i = 18r+ll, 1,..., 36r+ 

21}. Therefore, p(AV L ) np(AV^) = 0 for any i / j, 0 < i, j < 2r + 3. 
Since we have considered all additivities corresponding to the painting that 
could give rise to merging of components, it follows that p(A'P t ) for t = 
0,1,..., 2r+3 are 2(r+2) connected components. Furthermore, the intervals 
tg’ up] for i = 0, 1 ,..., q— 1 are all directly covered by the prescribed partial 
painting. □ 


Proof of Lemma 5.2. On the one hand, the function values (tto, h \, 
can be expressed in terms of (sq, si, •.., s r +i), as follows. 


i 717, 


T6 i 
7T6i+l 
7T6i+2 
7T6i+3 
^"6i+4 
*T6i+5 
^9r+5-3i 

7T9r+4-3i 

7T9r+3-3i 


6 Ey .| Sj + so — 2si — Si + 2sj+i, 

* = 0 , 1,... 

,r; 

6E}=i s j + so — 2si + 2sj+i, 

* = 0 , 1,... 

,r; 

6 E}=i + s o — 2sj + 3sj+i, 

* = 0 , 1,... 

,r; 

6 E}=i s j + s o ~ 2si + 4sj_|_i, 

* = 0,1,... 

,r; 

6 Ej=i s j + s o ~ 2si + 4sj+i + Sj+ 2 , 

* = 0,1,... 

,r- 1 

6 E}=i s j + s o — 2si + 5sj+i + Sj+ 2 , 

* = 0,1,... 

, r - 1 

9E, r = i^-3E;=i^ 

+ so — 2si + 7s r +i — Sj+i, 

* = 0,1,... 

,*■; 

9 Ej=i s j ~ 3 Ej=i s j 

+ so — 2si + 7s r +i — 2sj+i, 

* = 0,1,... 

,*■; 

9 Ej=i s j ~ 3 Ej=i s j 

+ sq — 2si + 7s r +i — 2sj+i — Sj-|_ 2 , 

* = 0,1,... 

,r - 1 
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(These formulas are obtained by integrating the slope values and are easily 
verified by induction.) By the symmetry condition, for t = 0,1,..., 9r + 5, 

^"9r+6+4 T Tg r _|_ 5 _£ — 18 f ( j ] Sj + 3so 6si T 14s r _|_i — 1. 

By the invariance property, for t = 0,1,..., 18r + 22, 7 j t = 7r q -t- 

On the other hand, the slope values (so, si, • • ■, s r +i) can be expressed 
in terms of the function values 7 Tf. so = ir\ — -ko, s\ = 712 — 7Ti, and st = 
7T 6t _ 8 - 7T 6t _9 for t = 2,3,..., r + 1. □ 


Proof of Lemma 5.3 . We show the ordering of (so, si, ..., s r +i) by consid¬ 
ering the subadditivity of ir. Since 7Ti + 7Ti > tt 2, we have so + so > so + si, 
and hence so > si. For i = 0, l,...,r — 1, the subadditivity condition 
T6i+3 T 7T9r+3— 3 j '7T9i—I-6-I-3* implies that 


i r i 

(eE Sj+so-2si+4sj + i)+(9 s j T s o — 2 s i+7s r + 1 —2 Si+ 1 —Si+ 2 ) 

7=1 7=1 7=1 

r r i 

£( 18 E Sj+3so—6si+14s r+ i)— (9 E s i- 3 E Sj + so—2s i+7s r + 1 —s j_)_ i ) . 

7=1 7=1 7 = 1 

After simplification, we have Si+i > Sj_)_ 2 - □ 
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